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Abstract. New detailed stellar yields of several elemen- 
tal species are derived for massive stars in a wide range of 
masses (from 6 to 120 Mq) and metallicities (Z =0.0004, 
0.004, 0.008, 0.02, 0.05). Our calculations are based on the 
Padova evolutionary tracks and take into account recent 
results on stellar evolution, such as overshooting and qui- 
escent mass-loss, paying major attention to the effects of 
the initial chemical composition of the star. We finally in- 
clude modern results on explosive nucleosynthesis in SNae 
by Woosley & Weaver 1995. The issue of the chemical 
yields of Very Massive Objects (from 120 to 1000 Mq) is 
also addressed. 

Our grid of stellar yields for massive stars is comple- 
mentary to the results by Marigo et al. (1996, 1997) on 
the evolution and nucleosynthesis of low and intermediate 
mass stars, also based on the Padova evolutionary tracks. 
Altogether, they represent a complete set of stellar yields 
of unprecedented homogeneity and self-consistency. 

Our new stellar yields are inserted in a code for the 
chemical evolution of the Galactic disc with infall of pri- 
mordial gas, according to the formulation first suggested 
by Talbot & Arnett (1971, 1973, 1975) and Chiosi (1980). 
As a first application, the code is used to develop a model 
of the chemical evolution of the Solar Vicinity, with a 
detailed comparison to the available observational con- 
straints. 

Key words: Stars: evolution, nucleosynthesis - stars: 
mass loss - stars: supernovas - Galaxy: Solar Neighbour- 
hood - Galaxy: evolution - Galaxy: abundances 
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1. Introduction 

In modelling the chemical evolution of our own Galaxy as 
well as of external galaxies, we deal with three basic ingre- 
dients: (1) the Star Formation Rate (SFR), (2) the Initial 
Mass Function (IMF) and (3) the stellar ejecta of differ- 
ent elements. Indeed, the chemical enrichment of galaxies 
is basically due to stars: while primordial nucleosynthesis 
produced mainly light elements, say up to 7 Li (Adouze 
1986, Walker et al. 1991), heavier species (metals) are due 
to stellar nucleosynthesis, either hydrostatic or explosive 
as in the case of supernovae (SNae). Therefore, a chemical 
model should estimate accurately the role of stars of dif- 
ferent masses in enriching the interstellar medium (ISM) 
with various elemental species. To this aim, one needs to 
consider the main nucleosynthetic processes taking place 
in stellar interiors, combined with a modelling of stellar 
evolution. 

New extended sets of stellar tracks for fine grids of 
masses and metallicities have recently become available 
from various research groups. The body of these mod- 
els shows that the initial chemical composition of a star 
bears on the details of its evolution. Therefore, we ex- 
pect that stellar ejecta depend not only on the mass of a 
star, but also on its metallicity. In spite of this evidence, 
most chemical evolution models rest upon grids of stellar 
yields derived for stars with solar composition, and only 
few examples of metallicity-dependent yields can be found 
in literature (Maeder 1992, hereinafter M92; Maeder 1993; 
Woosley & Weaver 1995, hereinafter WW95). 

In order to improve upon this, we derive here a wide 
and up-to-date set of ejecta in which the effect of the 
initial chemical composition is taken into account. We 
adopt the stellar models from the Padua library (Bressan 
et al. 1993; Fagotto et al. 1994a, b), which span the mass 
range from 0.6 to 120 Mq. Full sets are available for the 
metallicities Z=0.0004, Z=0.004, Z=0.008, Z=Z Q =0.02 
and .2=0.05, and we derive five corresponding sets of stel- 
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lar ejecta. Thus we get coherent prescriptions for stellar 
yields deduced from a unique and homogeneous grid of 
stellar models covering all mass ranges, instead of patch- 
ing together data from different sources. In addition, in 
the range of massive stars we take the effects of quiescent 
mass-loss starting from the ZAMS into account. We also 
point out that, unlike most stellar yields adopted in chem- 
ical models up to now, here we refer to evolutionary tracks 
calculated with convective overshooting. 

As a first application of our grid of stellar yields, we 
develop a chemical evolution model for the Solar Neigh- 
bourhood, which is the first template environment where 
nucleosynthetic yields and chemical models are to be cal- 
ibrated and tested. 

This paper is organized as follows. In Sec. 2 we briefly 
summarize the various groups in which stars of different 
mass can be classified according to the dominant physi- 
cal process and type of evolution. In Sec. 3 we deal with 
the role of mass loss in the evolution and nucleosynthe- 
sis of massive stars, and calculate the ejecta of the stellar 
winds predicted by our stellar models. In Sec. 4 we dis- 
cuss the effects of mass loss on the outcoming SNas and 
outline a recipe to calculate the explosive ejecta of SNae 
originating from our pre-SN models. Sec. 5 summarizes 
and comments upon the total ejecta of mass losing mas- 
sive stars with different initial metallicity. Sec. 6 describes 
the adopted prescriptions for the yields of low and inter- 
mediate mass stars. By means of simple semi-analytical 
prescriptions, in Sec. 7 we extend our grid of ejecta to 
Very Massive Objects up to 1000 M© (if they can form 
now or have existed in the past), since in some particu- 
lar astrophysical problems it can be useful to consider the 
possibility of chemical enrichment by these stars. Sec. 8 
introduces the mathematical formulation of our chemical 
model. In Sec. 9 we apply our model to the Solar Neigh- 
bourhood, discuss the space of parameters with respect to 
the available observational constraint and comment upon 
the chemical abundances and abundance ratios our model 
predicts. Finally, in Sec. 10 we draw some concluding re- 
marks. 

2. Definition of the mass intervals 

For the sake of better understanding, we summarize here 
how stars are classified according to the dominant phys- 
ical processes governing their evolution and fate. This is 
especially needed since we adopt models with convective 
overshoot for which the mass grouping is different as com- 
pared to standard models (cf. Chiosi 1986; Woosley 1986; 
Chiosi et al. 1992; and references therein). 

Low mass stars: They develop degenerate He-cores, un- 
dergo core He-flash, proceed to the formation of a 
degenerate CO-core, and after the AGB phase ter- 
minated by the loss of the envelope by stellar wind, 
become white dwarfs. In models with overshoot the 



upper mass limit (Mhcf) for this to occur is about 
1.6 -7- 1.7 Mq, depending on the chemical composition. 

Intermediate mass stars: They avoid core He-flash (he- 
lium is ignited in a non degenerate core), proceed to 
the formation of a degenerate CO-core, undergo the 
TP-AGB phase, and losing mass by stellar wind end 
up as white dwarfs. In models with convective over- 
shoot the upper mass limit (M up ) of this group is about 
5 Mq depending on the chemical composition. 

Quasi-massive stars: In models with convective overshoot, 
stars in the range 6 to 8 Mq undergo core C-burning in 
non degenerate conditions, but develop highly degener- 
ate O-Ne-Mg-cores. They become dynamically unsta- 
ble to e-capture instability and explode as supernovae. 

Massive stars: In the mass range 9 to 120 Mq there is 
the simultaneous occurrence of two phenomena: the 
dominant mass loss by stellar wind during the whole 
evolutionary history, and the completion of the nu- 
clear sequence down to the formation of an iron core 
in presence of strong neutrino cooling. The supernova 
explosion is triggered by core collapse induced by e- 
captures on heavy nuclei, photo-dissociation of Fe into 
a-particles, and rapid neutronization of the collapsing 
material. This is the range of classical Type II SNae. 
Either a neutron star of about 1.4 Mq or a black hole 
of larger mass is left over, depending on the efficiency 
of neutrino cooling during the previous stages. 

Very massive stars: Starting from the very initial stages, 
objects more massive than, say, 120 Mq are strongly 
pulsationally unstable, suffer from violent mass loss, 
and undergo pair- creation instability during core O- 
burning. The final outcome is regulated by the mass 
of the CO-core: at decreasing core mass, these stars can 
either collapse later to a black hole, or suffer complete 
thermonuclear explosion, or recover the behaviour of 
the massive stars above. 

3. Ejecta from the stellar wind of massive stars 

Massive stars contribute to the chemical enrichment of the 
ISM mainly by means of their final explosion in supernova 
(SN). Anyway, mass loss is also known to play a major role 
in the evolution of massive stars and stellar winds can 
have a non-negligible influence on stellar ejecta (Chiosi & 
Maeder 1986, hereinafter CM86; Chiosi & Caimmi 1979; 
Chiosi 1979; Maeder 1981, 1983, 1992, 1993). Mass loss 
has got two main effects on final ejecta: 

(1) stellar layers peeled off by stellar winds directly mix 
with and enrich the surrounding ISM, and evade further 
burning stages (direct effect); 

(2) mass loss affects the final mass Mf in of the star 
and the final masses of its He- and CO-core (M# e and 
Mco, respectively), thus affecting the resulting SN (indi- 
rect effect). 

In our adopted models, stars more massive than M > 
12 Mq are affected by stellar winds during their core H- 
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Table 1. Initial masses and corresponding final masses 



Fig. 1. Final mass vs. initial mass relation (in Mq) for massive 
stars of five metallicities (legend on top left). The dashed line 
corresponds to the case of constant mass 



and He-burning phases. All remaining evolutionary stages 
are so short-lived that mass loss has in practice no effect 
at all and the mass of star is frozen to the value with 
which it will undergo the final supernova explosion. This 
value of the mass is thereinafter referred to as the final 
mass Mfi n . The mass loss rates are included according to 
de Jager et al. (1988) from the ZAMS up to the de Jager 
limit; beyond such limit the star becomes a LBV candi- 
date and its mass loss rate is increased up to 10~ 3 M©/yr. 
When the surface hydrogen abundance falls below 0.3, 
the star is considered to enter its Wolf-Rayet stage and 
from then on the assumed mass loss rate follows Langer 
(1989). Since metallicity affects the bound-free and line 
opacity in the surface layers of massive stars, the effi- 
ciency of radiation pressure driven winds is expected to 
depend on metallicity according to M oc Z a , a ^0.5 
1.0 (Kudritzki et al. 1989; Kudritzki 1997; Leitherer & 
Langer 1991; Lamers & Cassinelli 1996). Evidence for the 
dependence of mass loss rates on metallicity comes from 
statistics of blue and red supergiants and of WR stars in 
our own Galaxy, in the Magellanic Clouds and in other 
galaxies of the Local Group (Maeder & Conti 1994). In 
our models, mass loss rates are scaled with metallicity ac- 
cording to Kudritzki et al. (1989, M oc Z - 5 ), therefore 
wind ejecta and the structure of SN precursors sensitively 
depend on the initial metallicity of the star, as shown by 
M92. For more details on the evolution of massive stars 
in presence of mass loss the reader is referred to CM86, 
Chiosi et al. (1992), Chiosi (1997). 
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In Fig. [l] we plot the final mass Mfi„ versus the ini- 
tial mass M for our five sets of models between 12 and 
120 M Q ; the corresponding values are listed in Tab. Q. 
As expected, higher metallicities correspond to lower fi- 
nal masses, due to larger mass loss rates. In particular, 
the most massive stars (M > 40 Mq) in the higher metal- 
licity sets end up with quite low final masses, typically 
4 v 6 Mq. These stars go through a WR phase, during 
which mass loss rates are high and strongly dependent on 
the stellar mass (M oc M 2 5 , Langer 1989); in this stage 
the stellar mass is decreasing rapidly, and the mass loss 
rate is also correspondingly decreasing, until all these stars 
reach very similar final masses. The final mass is very low 
especially in the case of high metallicity, because a more 
efficient wind during the Main Sequence phase allows the 
WR stage to be reached sooner and the intense WR mass 
loss to be active for a longer time. Indeed, to reach the 
final convergence mass the star needs to lose its H-rich en- 
velope soon enough (before the helium abundance in the 
core gets to Y ~ 0.5, according to Woosley et al. 1995). 

The ejecta due to stellar winds are calculated as 



j^wind 



where E wmd ; 



M(M,t)Xf(t)dt 



(1) 



iM is the amount (in Mq) of species i expelled in 
the wind by a star of mass M, M is the mass loss rate and 
Xf is the surface abundance of species i. The integration 
runs over the whole lifetime tm of a star of mass M. In 
practice, tm is the time when carbon ignites in the core, 
i.e. when the track ends. Wind ejecta are calculated for H, 
3 He, 4 He, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 0, 20 Ne, 22 Ne 
and are listed in Tab. 2. 

The amount and the chemical composition of the ex- 
pelled wind change with the mass and with the metallicity 
of the star; in order to explain this effect better, we con- 
sider stellar yields rather than ejecta. The stellar yield of 
a species i, pmi, is defined as the mass fraction of a star 
of mass M that has been newly synthesized as species i 
and then ejected (Tinsley 1980). If i has been converted 
into other species rather than produced, its yield piM is 
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Table 2. Ejecta (in M ) of the stellar wind for our model massive stars 
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negative. On the basis of this definition, the total amount 
of species i ejected by a star of mass M is: 

E iM = (M - M r ) X? + Mp iM (2) 

where M r is the remnant mass and X® is the initial abun- 
dance of species i in the star when it was born. The first 
term is the amount of "original" species i eventually ex- 
pelled, while the second term is the positive or negative 
contribution of newly synthesized i. Therefore, yields allow 
us to estimate the nucleosynthetic production due to the 
star itself, by distinguishing the newly synthesized frac- 
tion from the fraction that was already present in the gas 
whence the star formed. The stellar yields of the wind for 



our set of masses and metallicities are calculated according 
to: 

i-tm 

Mp^ d = / M(M, t) [Xf{t) - Xf\ dt (3) 
Jo 

The yields formalism is basically due to Tinsley (1980), 

but the previous formulae are taken as revised by M92. 

The wind yields we obtain are listed in Tab. 3 and can be 

commented on as follows (see also M92). 

1 H: Hydrogen yields are always negative and increase in 
modulus with increasing mass and metallicity, since 
the more efficient mass loss is, the larger amount of 
new helium and metals (synthesized at the expense of 
1 H) is revealed on the surface and released in the wind. 

3 He: The abundance of 3 He increases at the very outer 
edge of the H-burning region, due to the first reactions 
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Table 3. Yields pTm 111 (see Eq. 3) of the stellar wind for our massive stars. When the amount of a newly synthesized/consumed 
element is less than 1% of its initial one, the corresponding yield is assumed to be not significative, i.e. is set equal 
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of the p-p chain; then 3 He is rapidly turned to 4 He. Its 
yields are generally negative, except in some cases of 
low mass loss where the 3 He-enriched layers are lost, 
rather than inner ones. 

4 He: Helium yields are larger for larger masses and, at 
any given mass, for higher metallicities, because strong 
winds can take large amounts of new 4 He away. 

12 C — le O: Carbon yields are negative in most cases, since 
the CNO cycle turns 12 C into 14 N. Only when mass loss 
gets very efficient, i.e. for high masses and metallici- 
ties, carbon yields can be positive because He-burning 
products can be revealed on the surface (WC and WO 
stars). Quite the same holds for oxygen, whose yields 
are seldom positive. If mass loss is extreme, anyway, 
12 C and 16 O yields may even decrease because most of 
the mass is rapidly lost in the wind in the form of 4 Hc. 

14 N — 13 C — iy O: In our models of massive stars, nitrogen 
is produced only as a secondary element, therefore its 
yields sensitively grow with metallicity. The same holds 
for 13 C and for 17 0, which are secondary products of 
the CNO cycle as well. 

15 N: 15 N is quickly destroyed in the CNO cycle, therefore 
its yields are always negative (although small) and in- 



creasing in modulus with metallicity, as the efficiency 
of the CNO cycle increases. 

18q. isq j g destroyed in the CNO cycle; it is later pro- 
duced by a capture on 14 N during He-burning, but 
then it is turned into 22 Ne by a new a capture, so 
that the abundance peak of 18 is very thin (CM86). 
Therefore, even when mass loss is so efficient as to re- 
veal He-burning products on the surface, 18 O yields 
usually remain negative. 

22 Ne: 22 Ne yields are negligible for relatively low masses 
and/or metallicities, while when mass loss gets very 
efficient there is a sensitive ejection of new 22 Ne as a 
by-product of He-burning remained uncovered on the 
surface of WR stars. 

20 Ne: 20 Ne yields are always negligible because, although 
20 Ne can be produced by a capture on 16 O during He- 
burning, this holds for so advanced stages that even an 
extreme mass loss can't peel off the layers where 20 Ne 
is synthesized. 

4. The supernova explosion 

To determine the contribution of quasi-massive and mas- 
sive stars in the chemical enrichment of the ISM wc need 



(i 
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Table 4. He-core mass (in Mq) of model stars in four different evolutionary stages (see text). Here Mhb is defined as the 
mass-point where the hydrogen content falls to 
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to know the ejecta of the final SN explosion, which ac- 
tually plays the major role. Although mass loss is well 
known to have substantial effects in the evolution of mas- 
sive stars, SN models usually start from constant mass pre- 
SN structures. Only few examples are found in literature 
of stellar models with mass loss followed up to the very 
final stages and the explosion (e.g. Woosley et al. 1993, 
1995), but no extended studies exist yet of explosive nu- 
cleosynthesis of mass losing massive stars for wide sets of 
masses and metallicities. To get the global nucleosynthetic 
production for these stars one is generally forced to link, 
somehow, pre-SN models evolved including stellar winds 
and SN models based on constant mass calculations. Since 
our tracks stop at the beginning of C-ignition in the core, 
here we need to perform such a link as well. 

Under the hypothesis that the stellar core, which drives 
the explosion mechanism, evolves uncoupled from the en- 
velope, SN models were calculated letting He-cores, i.e. 
pure helium stars, evolve (Arnett 1978, 1991; Nomoto 
& Hashimoto 1986; Thielemann et al. 1996; Woosley & 
Weaver 1986). These bare helium stars are later identi- 
fied with the helium core (Mjie) of stars with initial, total 
mass M. By means of a suitable relation between the to- 
tal mass M of the star and its Mff e , the ejecta obtained 
for bare He-cores were then shifted into a scale of stel- 
lar masses. The envelope layers, and also the wind con- 
tribution in the case of mass losing stellar models, were 
added to the ejecta of the corresponding He-core (Chiosi 
& Caimmi 1979; Chiosi 1979; Baraffe & El Eid 1991). Yet, 
as pointed out by Maeder (1981, 1984, 1992), this method 
leads to some inconsistency, because a stellar core does 
not necessarily evolve like an isolated one. Unlike a bare 
He-core, its size may change during stellar evolution: it 



can grow or recede in time, depending also on dredge- 
up episodes and mass loss. Table f| shows this effect by 
displaying Mjje values in four evolutionary stages of our 
tracks: (1) when H is exhausted in the centre, (2) when He 
ignites in the core, (3) when He is exhausted in the centre 
and (4) when C ignites in the core. Mn e gets growing in 
the course of evolution for less massive and/or low metal- 
licity stars; in stars where mass loss is very efficient, Mff e 
falls in later stages. 

Therefore, it's hard to find an univocal relation be- 
tween the total initial mass of the star and M# e . The 
meaningful value of Mn e for the sake of the resulting SN 
might be the one corresponding to the end of core He- 
burning or to C-ignition, since following stages are so fast 
that mass loss cannot alter the structure of the star any 
more. However, even the "final" He-core mass may be in- 
adequate to fix the resulting explosion. In particular, in 
the WR stage the strong mass loss proceeding after the 
He-core is uncovered can alter the evolution of the core 
and the final nucleosynthesis. Woosley et al. (1993) fol- 
low the evolution of a mass losing 60 Mq star evolving 
through a WR stage and ending as a 4.25 Mq bare He- 
core at the time of explosion. Such a structure is compared 
to a 4.25 M & pure He star evolved without mass loss. The 
final thermal and dynamical structure of the two objects 
results very similar, but their chemical structure is differ- 
ent: although the mass of the pre-SN He-core is the same, 
Mco is different and this affects all subsequent burning 
stages, the resulting entropy distribution and the final ex- 
plosion. 

Therefore, here we prefer to choose Mco as the char- 
acterizing parameter for the resulting SN, as suggested by 
M92. Indeed, C-burning and following burning stages are 
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Fig. 2. He-core mass vs. initial mass relation (in M ) for 
massive stars of five metallicities. Asterisks indicate He-core 
masses from Woosley & Weaver 1986, relevant to constant mass 
pre-SN models 



Fig. 3. CO-core mass vs. initial mass relation (in M ) for 
massive stars of five metallicities. Asterisks indicate CO-core 
masses from Woosley & Weaver 1995, relevant to constant mass 
pre-SN models 



so quick that mass loss hasn't got time enough to influ- 
ence inner cores and can be ignored (see also Woosley 
et al. 1993, 1995). Mco at carbon ignition well corre- 
sponds to Mco at the time of explosion. Basing on Mco 
we perform a link with SN models for all our massive 
stars, i.e. for our models with M > 6 M . Even though our 
models include quiescent mass loss only for M > 12 M Q , a 
link with SN models is made also for 6 M Q < M < 12 M , 
since in our models with overshooting M up ~ 5 M and 



all higher masses are expected to undergo SN explosion. In 
general, we remind that due to overshooting any given M 
here corresponds to larger M# e and Mco with respect to 
"standard" models, and since the final fate of a star is ul- 
timately related to the size of its core the significant mass 
ranges implying different kinds of SN arc shifted down- 
ward in mass. 



Table 5. Initial masses and corresponding He-core masses Table 6. Initial masses and corresponding CO-core masses 
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25.40 


17.55 


10.65 


5.77 


3.98 


100 


49.16 


45.47 


13.83 


6.96 


3.92 


120 


56.67 


38.43 


13.95 


6.44 


3.93 



M 


Mco 




Z=0.0004 


Z=0.004 


Z=0.008 


Z=0.02 


Z=0.05 


6 


1.11 


1.11 


1.05 


1.01 


1.00 


7 


1.28 


1.20 


1.22 


1.10 


1.16 


9 


1.80 


1.83 


1.75 


1.68 


1.71 


12 


2.78 


2.82 


2.73 


2.64 


2.79 


15 


3.97 


3.96 


3.89 


3.89 


4.19 


20 


6.25 


6.15 


6.02 


5.94 


6.84 


30 


11.23 


10.95 


10.86 


11.02 




40 


16.39 


16.03 


15.70 


4.59 


2.95 


60 


24.58 


16.40 


9.97 


5.18 


3.41 


100 


48.45 


44.94 


13.29 


6.39 


3.36 


120 


49.51 


35.28 


13.40 


5.87 


3.35 
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Table 7. Total ejecta (in Mq) for quasi- massive stars of 6 and 7 Mq, for different metallicities 



Z— 0.0004 



M 


*H 




4 Hc 


12 c 




i4 N 


N 


16 G 




18 Q 






M r 


6 
7 


3.22 
3.68 


1.70E-04 
1.45E-04 


1.48 
2.00 


2.42E-04 
2.67E-04 


1.48E-05 
1.71E-05 


5.23E-04 
7.25E-04 


1.82E-07 
2.02E-07 


8.12E-04 
9.10E-04 


3.94E-06 
4.37E-06 


1.35E-06 
1.49E-06 


1.93E-05 
2.34E-05 


5.13E-06 
5.99E-06 


1.30 
1.30 


Z=0.004 


M 


% 


;i Hc 


4 Hc 


12 C 


13 C 


14 N 


Ib N 


16 G 


A 'o 


18 


^Ne 


^Nc 


M r 


6 
7 


3.16 
3.61 


1.S0E-04 
1.64E-04 


1.52 
2.07 


2.80E-03 
3.22E-03 


1.65E-04 
1.99E-04 


4. 72 E- 03 
6.42E-03 


2.04E-06 
2.34E-06 


8.17E-03 
9.35E-03 


3.99E-05 
3.93E-05 


1.56E-05 
1.80E-05 


1.93E-04 
2.34E-04 


4.73E-08 
5.45E-08 


1.30 
1.30 



M 






4 Hc 


12 C 


13 C 


14 N 


iD N 


16 




1H 






M r 


6 


3.15 


1.97E-04 


1.51 


5.97E-03 


3.30E-04 


8.53E-03 


4.43E-06 


1.69E-02 


1.09E-04 


3.33E-05 


3.86E-04 


1.03E-04 


1.30 


7 


3.59 


1.80E-04 


2.06 


6.73E-03 


3.98E-04 


1.21E-02 


4.93E-06 


1.91E-02 


1.25E-04 


3.79E-05 


4.68E-04 


1.20E-04 


1.30 



M 


*H 


3 He 


4 Hc 


12 c 


13 c 


14 N 




16 G 


Wo 


ia 






M r 


6 
7 


3.02 
3.44 


2.46E-04 
2.34E-04 


1.58 
2.15 


1.60E-02 
1.79E-02 


8.49E-04 
1.00E-03 


1.84E-02 
2.72E-02 


1.18E-05 
1.32E-05 


4.38E-02 
4.97E-02 


5.24E-04 
5.71E-04 


8.83E-05 
9.89E-05 


9.64E-04 
1.17E-03 


2.57E-04 
2.99E-04 


1.30 
1.30 


Z— 0.05 


M 


i H 


:i Hc 


4 Hc 


12 C 


13 C 


14 N 




W Q 


A 'o 


18 Q 


^Ne 


^Nc 


M r 


6 
7 


2.49 
2.88 


3.37E-04 
3.59E-04 


1.97 

2.53 


4.03E-02 
4.61E-02 


2.08E-03 
2.51E-03 


4.39E-02 
6.13E-02 


3.02E-05 
3.38E-05 


0.111 
0.129 


1.66E-03 
1.99E-03 


2.21E-04 
2.53E-04 


2.41E-03 
2.92E-03 


6.42E-04 
7.49E-04 


1.30 
1.30 



4--1- The adopted Mhc(M) and Mco(M) relations 

In Fig. H we plot the relation between the initial mass M 
and the corresponding He-core mass Mjj e at the end of 
the evolutionary track. For the sake of comparison, aster- 
isks indicate A/# e values given by Woosley & Weaver 1986 
for pre-SN models calculated at constant mass. Up to ~ 
30 Mq the He-cores are more massive in our models due 
to overshooting, while for larger masses the effects of mass 
loss dominate. Indeed, in the range 10 4- 30 Mq stellar 
evolution is more sensitive to the treatment of convection, 
while for larger masses the assumptions about mass loss 
have an overwhelming influence (CM86). 

In Fig. H we plot Mco at the end of the track versus 
M; in our models, we define Mco as the mass interior to 
the outer edge of the He-burning region. Asterisks indicate 
Mco values deduced from the recent SN models by WW95 
(see App. A), calculated from pre-SN structures with no 
mass loss and with no overshooting. 

In Tab. |^ and Tab. || we list the values of M# e and 
of Mco respectively, for our quasi-massive and massive 
stars (M > M up —5 Mq). In models with overshooting 
all these stars eventually give SNae, and their core masses 
are prerequisite to characterize the outcoming supernova 
and its ejecta. 



4-2. Electron capture supernovas 

Quasi-massive stars (8 -j- 10 Mq in standard models and 
6-^8 A/© in models with overshoot) generally develop a 
degenerate O-Ne-Mg core after C-burning and eventually 
explode as "electron capture" supernovas (EC SNae), leav- 
ing a neutron star of ~1.3 Mq as a remnant and expelling 
very few heavy elements (Nomoto 1984, 1987; Mayle & 
Wilson 1988). Recent calculations seem to confirm that 
a degenerate O-Ne-Mg core reaching 1.375 Mq, the limit 
mass for electron captures on 24 Mg and 20 Ne, develops 
density conditions extreme enough that it collapses to a 
neutron star, rather than undergo a thermonuclear explo- 



sion due to disruptive O-Nc burning in degenerate mate- 
rial (Gutierrez et al. 1996). 

But the final stages of stars in this mass range are still 
debated. These stars develop a double-shell structure sim- 
ilar to that of lower mass stars, and experience a thermally 
pulsing, mass losing phase as well. Nomoto's models didn't 
include thermal pulses nor mass loss, but the detailed evo- 
lution of a model 10 Mq star has been recently followed by 
Garcia-Berro & Iben (1994), Ritossa et al. (1996). Their 
models suggest that such a star might be unable to reach 
the limit core mass for electron capture and explode in 
an EC SN: it might rather get rid of the overlying layers 
through a superwind phase typical of TP-AGB stars and 
result in a massive (1.26 Mq) O-Ne-Mg white dwarf. 

By comparing the size of our stellar cores with those 
of Nomoto (1984, 1987) and Ritossa et al. (1996), in our 
models with convective overshooting such a behaviour is 
expected between 6 and 8 Mq (M C o ~ 1.0 -i- 1.3 A/ 
at the beginning of core C-burning). For our models in 
this mass range, we assume that the inner 1.3 Mq remain 
locked in the remnant, either a neutron star or a white 
dwarf, while the overlying layers are expelled, either by 
a final explosion or by a superwind during the TP-AGB 
phase. The resulting ejecta for this mass range are pre- 
sented in Tab. (?]. 

We thus neglect the products of explosive nucleosyn- 
thesis in case of SN explosion, which is reasonable since 
EC SNae are believed to produce negligible amounts of 
heavy elements (not more than 0.002 Mq, Mayle & Wil- 
son 1988). More important, we are neglecting the effects of 
thermal pulses and of the III dredge-up in the nucleosyn- 
thetic yields of SAGB stars (SAGB = Super- Asymptotic 
Giant Branch, see Garcia-Berro & Iben 1994, Ritossa 
et al. 1996). No full calculations of stellar yields for this 
kind of stars exist yet in literature, neither for standard 
models nor for models with overshooting. Further investi- 
gation of nucleosynthesis and mass loss in this mass range 
should be done. 
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Table 8. Remnant mass M r and ejecta E?° (in M ) of CO -cores for different elements, as deduced from WW95 



Z=Z B (case S in WW95) 



Mco 


M r 


L'2 C 


ib N 


1B 


a) Nc 




•2H si 


■V2 S 


4U Ca 




1.49 


1.26 


3.13E-02 


1.29E-05 


6.47E-02 


1.58E-02 


4.34E-03 


1.55E-02 


5.94E-03 


7.50E-04 


6.93E-02 


1.72 


1.30 


5.59E-02 


3.77E-05 


0.130 


6.55E-03 


2.92E-03 


8.42E-02 


7.16E-02 


1.38E-02 


4.34E-02 


1.95 


1.40 


8.95E-02 


4.01E-06 


0.192 


2.67E-02 


1.07E-02 


5.13E-02 


2.16E-02 


2.95E-03 


0.133 


2.53 


1.40 


0.134 


9.20E-05 


0.591 


9.08E-02 


2.03E-02 


0.102 


5.85E-02 


1.03E-02 


0.115 


3.54 


1.70 


0.215 


9.34E-05 


1.02 


0.254 


4.78E-02 


0.128 


4.89E-02 


5.95E-03 


6.59E-02 


4.16 


1.90 


0.250 


6.83E-05 


1.32 


8.20E-02 


1.75E-02 


0.267 


0.123 


1.26E-02 


0.100 


4.73 


2.00 


0.178 


8.48E-05 


1.83 


8.03E-02 


2.34E-02 


0.278 


0.146 


1.34E-02 


8.84E-02 


5.37 


1.90 


0.203 


1.10E-04 


2.26 


4.33E-02 


3.30E-02 


0.345 


0.165 


1.63E-02 


0.205 


6.59 


1.90 


0.281 


1.64E-04 


3.12 


3.64 


9.65E-02 


0.303 


0.134 


1.56E-02 


0.129 


8.67 


4.10 


0.241 


9.62E-05 


3.51 


0.346 


0.201 


8.38E-02 


5.35E-03 


3.25E-05 


4.80E-05 


11.12 


7.20 


0.253 


5.13E-05 


2.92 


0.476 


0.111 


2.60E-02 


1.34E-03 


2.97E-05 





15.41 


12.30 


0.267 


2.07E-05 


2.21 


0.408 


5.53E-02 


8.64E-03 


4.62E-04 


6.95E-06 













Z=Q.l Z (casc I 


' in WW95) 










Mco 


M r 


VJ C 


i0 N 


lb Q 


2u No 




' 2 *Si 




40 Ca 


bb Fo 


1.81 


1.30 


8.66E-02 


9.74E-06 


0.137 


1.04E-02 


4.95E-03 


2.95E-02 


1.32E-02 


2.78E-03 


0.177 


2.01 


1.20 


0.106 


1.65E-05 


0.282 


4.04E-02 


1.63E-02 


6.43E-02 


3.11E-02 


5.16E-03 


0.180 


2.52 


1.40 


0.139 


4.79E-05 


0.546 


7.26E-02 


2.11E-02 


7.03E-02 


3.35E-02 


5.83E-03 


0.194 


3.43 


1.60 


0.222 


5.82E-05 


0.983 


0.204 


7.71E-02 


0.126 


4.74E-02 


7.75E-03 


0.136 


4.32 


1.90 


0.242 


7.17E-05 


1.51 


9.36E-02 


5.07E-02 


0.231 


0.128 


1.86E-02 


0.121 


5.11 


2.00 


0.228 


1.01E-04 


2.11 


3.32E-02 


1.64E-02 


0.325 


0.183 


2.27E-02 


0.123 


6.06 


1.90 


0.263 


1.25E-04 


2.89 


4.13E-02 


2.95E-02 


0.394 


0.216 


2.16E-02 


0.202 


8.06 


2.60 


0.299 


1.28E-04 


4.10 


0.634 


0.319 


7.06E-02 


1.95E-03 


5.06E-05 





10.55 


6.50 


0.314 


5.89E-05 


3.08 


0.488 


0.144 


2.28E-02 


8.32E-04 


5.06E-05 





12.65 


8.90 


0.332 


3.55E-05 


2.70 


5.99 


7.61E-02 


6.02E-03 


5.27E-04 


5.51E-05 













Z=0.01 Z (case 


r in WW95) 








Mco 


Mr 


V2 C 


io N 


lb 


aJ Nc 




28 si 


32 g 


4u Ca 


bb Fo 


1.80 


1.30 


8.98E-02 


6.65E-06 


0.141 


9.46E-03 


5.75E-03 


3.12E-02 


1.40E-02 


2.61E-03 


0.152 


2.01 


1.40 


0.109 


1.01E-05 


0.209 


2.48E-02 


1.05E-02 


4.14E-02 


2.18E-02 


4.89E-03 


0.191 


2.42 


1.50 


0.149 


2.14E-05 


0.422 


4.00E-02 


2.12E-02 


7.32E-02 


3.46E-02 


6.46E-03 


0.171 


3.18 


1.50 


0.194 


4.76E-05 


0.951 


0.185 


3.58E-02 


0.101 


5.60E-02 


9.94E-03 


0.135 


4.31 


1.90 


0.208 


8.71E-05 


1.62 


5.46E-02 


1.66E-02 


0.244 


0.148 


1.84E-02 


7.30E-02 


4.83 


1.90 


0.248 


9.57E-05 


1.94 


4.11E-02 


5.73E-02 


0.311 


0.148 


1.93E-02 


9.88E-02 


5.84 


1.70 


0.279 


1.27E-04 


2.83 


4.37E-02 


3.95E-02 


0.385 


0.205 


1.94E-02 


0.187 


8.10 


3.10 


0.315 


1.35E-04 


3.98 


0.357 


0.255 


0.116 


3.24E-03 








10.20 


5.20 


0.343 


8.29E-05 


3.67 


0.711 


0.180 


1.78E-02 


1.39E-04 








12.39 


8.90 


0.340 


3.82E-05 


2.61 


0.532 


7.41E-02 


3.60E-03 












4-3. Iron-core collapse supernova 

In the range M >8 Mq, our models go through the whole 
sequence of nuclear burning stages and eventually explode 
after iron-core collapse. Very few models explode as "pair 
creation" SNae (PC SNae) before building an iron core (see 



4.4). To get the final outcome of iron core collapse SNae, 



we link our models to the recent SN models by WW95, 
which cover a wide range of masses and metallicities (from 
11 to 40 M Q and from Z=0 to Z=Z Q =0.02). In WW95, 
pre-SN models are calculated at constant mass and con- 
vective regions are treated adopting the Schwarzschild cri- 
terion and semiconvection; the models are exploded by 
letting a piston move outward, tuning it so that the fi- 
nal ejecta gain a typical kinetic energy at infinity (1.2 x 
10 51 erg, although for M > 30 Mq higher energies are also 
considered in cases B and C). 

WW95 give the total ejecta, i.e. the ejecta relevant to 
the whole star. In order to link the outcome of these SNae 
to our models, we need to discriminate the contribution of 
the CO core alone in the models of WW95. We consider 
the following elemental species: H, 4 He, 12 C, 13 C, 14 N, 
15 N, 16 0, 17 0, 18 0, 20 Ne, 22 Ne, 24 Mg, 28 Si, 32 S, 40 Ca, 
56 Fe. By means of the method described in App. A, for 
various metallicities referred to in WW95 we derive suit- 



able relations between the mass of the CO-core and the 
amount of different elements expelled by the core in the ex- 
plosion. These Ef° — Mco relations are shown in Tab. || 
and in Fig. ||. We notice quite a regular behaviour of the 
remnant mass and of the ejecta of most elements as a func- 
tion of Mco > even with varying initial metallicity. This 
confirms that Mco is a good parameter for our link (but 
see also App. B). 

By means of a linear interpolation with respect to Mco 
and to metallicity in such relations, we get the ejecta Ef° 
from the CO-cores of our models. Then we add (1) the 
composition of the overlying layers still left in our models 
at the time of explosion, and (2) the ejecta E™ md from the 
stellar wind (Tab. 4), to get our total ejecta as the sum: 



Ei = E 



CO 



E° 



E]' 



(4) 



4-4- Pair creation supernova 



Massive stars with larger cores (Mjie ^ 35 Mq) are 
known to undergo pair creation events during O-burning. 
This "pair creation" instability (Fowler & Hoyle 1964) 
may lead either to violent pulsational instability with ejec- 
tion of some external layers and later iron core collapse 
(35 Mq < MHe % 60 Mq), or to complete thermonuclear 



10 
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Fig. 4. Remnant mass and ejecta of the CO-core as functions of Mco for three different metallicities, as deduced from WW95 
(see text). Circles: case Z=Zq; squares: case Z—Q.l Zq \ triangles: case Z=Q.Q1 Z Q . 
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explosion (60 Mq < Mjje < HO Mq), or to collapse to 
black hole after He-exhaustion (Mf/ e ^ HO Mq); for de- 
tails refer to Woosley (1986). In terms of CO-core masses, 
up to Mqo^ 57 Mq we expect pulsational instability, in 
the range Mco=57 -j- 106 M a complete disruption is 
supposed to occur, while larger cores collapse to black 
holes. 

Because of mass loss, in our stellar models we seldom 
face such large cores and "pair creation" supernova? (PC 
SNae) are likely to occur only in the case of the higher 
mass, low metallicity stars (M=100-120 M Q , Z=0.0004- 
0.004), as mass loss in not very efficient. These stars 
have He-cores in the range 35 -j- 60 M (Tab. ; follow- 
ing Woosley (1986) we assume that they expel some of 
the external layers when pair instability develops, while 
whatever remains bound finally collapses to a black hole. 
In order to establish what is to be expelled in our CO- 
cores and what is to be locked in the remnant, we in- 
terpolate between the cases Mjy e =45 Mq and 55 M Q 
of Woosley (1986), corresponding to Mco ~ 41 M and 
Mco ~ 52 Mq, respectively. For Mco=41 Mq only over- 
lying layers are expelled, while for the Mco =52 Mq the 
inner 27 Mq remain bound. Detailed ejecta of this case, 
as deduced from Tab. 18 of Woosley (1986), are shown in 
Tab. [I The total ejecta of our stellar models experiencing 
a PC SN event are calculated according to Eq. ([|). 
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Fig. 5. Fractional remnant mass KM=M r /M and total 
yields piM for some elements vs. the initial stellar mass M, 
for the Z=0.0004 set. Solid line: RM; dashed line: helium 
yields; dotted line: carbon yields; dash-dotted line: oxygen 
yields;long-dashed line: heavy elements yields. 



Table 9. I s * column: M He values of PC SNae models of 
Woosley (1986). 2 nd column: corresponding CO-core masses. 
3 r to 6 th column: ejecta from the CO-core for significant el- 
emental species. 7 th column: remnant mass. All quantities are 
expressed in Mq 





Mco 


V2 C 


16 






M r 


45 


41 














41 


55 


52 


1.1 


21 


1.5 


1.1 


27 



5. The total ejecta of massive stars 

The total ejecta for our stellar models, obtained as out- 
lined in the previous sections, are listed in Tab. [l^. 

In Figures ^| to ^ we plot the fractional remnant mass 
M r /M and the total yields PiM for some representative 
elements vs. the stellar mass M, for our five sets of metal- 
licity. A trend with metallicity is evident. 

In the lower metallicity sets (Z=0.0004 and Z=0.004) 
remnant masses are very large due to lower mass-loss ef- 
ficiency: massive CO-cores are built in the pre-SN phase, 
resulting in high remnant masses (i.e. black holes, cfr. 
Tab. ||). We notice an anti-correlation between the rem- 
nant mass and the oxygen yield: for instance, in the case 
Z=0.0004, M=40 and 60 M© (M C o~ 18-24 M ) almost 
the whole CO-core collapses to a black hole and a very 



small amount of new oxygen is released. For larger masses 
PC SNae are found and oxygen yields increase again, while 
remnant masses decrease. This is in line with what Langer 
& Woosley (1996) suggest: at low metallicity (i.e. in the 
early phases of galactic evolution) stars in the range 30 -j- 
100 Mq may form black holes and have a secondary role 
in chemical enrichment, while higher masses may give an 
important contribution through PC SN33, also depending 
on the IMF in the early epochs. For solar metallicity, on 
the contrary, mass loss inhibits black hole formation and 
PC SNae events. 

In the case Z=0.008, the effects of mass loss are getting 
evident: for large masses (M > 40 Mq) remnant masses 
are sensibly smaller than in the previous sets, while helium 
and carbon yields arc larger because of the contribution 
from the wind. 

These effects are exalted in the higher metallicity sets 
(Z=Zq=0.02 and Z=0.05), where large initial masses end 
as WR stars with small masses and small cores (cfr. Tab. 
and Tab. ^), resulting in relatively low oxygen yields but 
very large helium yields. Carbon yields decrease again 
moving from the Z=0.008 to the Z=0.Q5 case, both be- 
cause a higher mass loss rate is able to take more helium 
away before it is turned to carbon, and because with an 
increasing metallicity an increasing fraction of the original 
carbon is turned to 14 N in the CNO cycle. 
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Table 10. Total ejecta (in M©) for model stars between 9 and 120 M© for different metallicities 



Z=0.0004 





1 H 




4 Hc 


12 c 


13 c 


14 N 


lb N 




lv o 


18 




— — 








!20tvt 




N 


e 


^ 4 M g 


2H" Si 


sa s 


4,) Ca 




M r 


9 


4.57 


1 


22E 


04 


2.63 


8.97E-02 


2 


27E 


05 


9 


31E-04 


7 


16E-06 




0.140 


4 


78E 


06 


2 


43E-05 
















9.43E-03 


2 


96E 


05 




61E-03 


3 


08E-02 


i 


38E-02 


2 


60E 


03 




0.155 


1.35 


12 


5.77 


9 


50E 


05 


3.45 


0. 170 


2 


91E 


05 


1 


28E-03 


3 


58E-05 




0.671 


5 


09E 


06 


2 


97E-05 
















0. 108 


3 


82E 


05 


2 


91E-02 


8 
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As far as heavy elements ( Si 



3 Ce 



Fe) are 



concerned, the bulk of contribution seems to come from 
stars between 10 and 20 Mq. For larger masses large rem- 
nants are found, so that most heavy elements produced 
remain locked in the remnant, at least for Z < 0.008. In 
the case of solar or super-solar metallicity, also stars with 
M > 20 M Q result in low CO-cores and low-mass rem- 
nants and can release some heavy elements. 

6. Low and intermediate mass stars 

For low and intermediate mass stars (M < 5 Mq) we 
adopt the yields by Marigo et al. (1996, 1997), which are 
calculated basing on the same series of stellar tracks we 
adopt for massive stars. Specifically, by means of a semi- 
analytical model Marigo et al. (1996) extend the tracks of 
the Padua library, which stop at the end of the E-AGB 
phase, to the TP-AGB and III dredge-up phase. Further 
improvements of the semi-analytical formulation also al- 
lows to properly take into account the effects of envelope 
burning on the yields of stars between 4 and 5 Mq (Marigo 
et al. 1997). The related production of primary 14 N has re- 
markable effects upon the predictions of chemical models 



for the Solar Neighbourhood (§ 9.7). 

Thus we extend our grid of stellar ejecta to the range 
of low and intermediate mass stars and provide a complete 
data-base for the chemical models of galaxies. This data- 
base is referred to a unique grid of stellar models for the 
whole range of stellar masses, which gives it coherence and 
homogeneity of basic physical prescriptions. 



Table 11. Lifetimes, H-burning core masses and CO-core 
masses for Very Massive Objects 
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7. Very Massive Objects 

The structure and evolution of very massive massive 
objects (VMOs), i.e. stars in the mass range 120< 
M <1000 Mq, was explored in past years, mainly in the 
framework of a primeval Population III of very massive 
stars formed by metal-free gas (Bond et al. 1984; Carr 
et al. 1984; El Eid et al. 1983; Ober et al. 1983; Woosley 
& Weaver 1982). Population III was invoked to solve the 
G-dwarf problem and to explain the non-zero metallicity 
of Population II stars, to provide a substantial amount of 
dark matter by means of black hole formation, to account 
for black holes in AGNs, to explain the reionization of the 
Universe, to produce primordial helium and to account for 
some difficulties in Big Bang nucleosynthesis; for a review 
see Bond (1984). 

Some recent interest toward stars more massive than 
the canonical limit of say 120 Mq come from studies of 
chemical evolution and population synthesis in galaxies 
under unconventional IMFs (Chiosi et al. 1997) such as 
that proposed by Padoan et al. (1997) which formally al- 
lows for a non-zero, though small, probability that objects 
more massive than 120 Mq can form. 

We are not going to discuss here the possible exis- 
tence of VMOs in present-day or primeval galaxies, nor 



the cosmological consequences of Population III; rather, 
we want to make reasonable assumptions about nucleosyn- 
thesis and mass loss in VMOs and to extend our five grids 
of ejecta up to 1000 Mq, for the purpose of future use 
if required by the particular problem under examination. 
Notice that we do not include chemical enrichment from 
VMOs in our chemical model of the Solar Neighbourhood 
(see Sect. 9). 

VMOs are pulsationally unstable during H-burning, so 
once they form they are likely to undergo violent mass 
loss until they fall below the critical mass for pulsational 
stability; this limit is generally around 100 Mq, though 
depending on metallicity. We assume that VMOs lose mass 
at a rate of 10~ 3 M Q /yr (comparable to the extreme mass 
loss rate adopted in our stellar tracks beyond the de Jager 
limit), independently of metallicity since mass loss here 
is driven by instability rather than by radiation pressure. 
We stop this phase of paroxysmal mass loss when the mass 
has decreased down to 150 Mq. 

Meanwhile, we assume that the star has been burning 
hydrogen in the inner 50% of its mass (He-cores of VMOs 
involve ^0.56 of the initial mass, see the references quoted 
above) and that the hydrogen content X in the core has 
been linearly decreasing with time, consistently with the 
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behaviour of our high-mass stellar tracks. The hydrogen 
content in the core at any time t is given by: 



X(t) = X° ( 1- 



t 

TH 



(5) 



where X° is the initial hydrogen abundance in the star and 
th is the H-burning timescale, assumed following the ref- 
erences above and listed in the second column of Tab. [ll[ 
At the end of the violent mass loss phase (i.e. when 
the current mass falls below 120 M Q ) we adopt a rate for 
radiation pressure driven mass loss derived from Lamers 
& Cassinelli (1996): 



M = 2.11 x 1CT 14 x M- 4A21 



(6) 



scaled with metallicity in conformity with our stellar 
tracks for massive stars (M oc Z - 5 ). This mass loss rate 
holds until the star becomes a WR star; in case it enters 
this stage, for the rest of its lifetime we adopt the mass 
loss rate by Langer (1989), scaled with metallicity. For 
the sake of consistency with our stellar tracks, we assume 
that the WR stage begins if the surface hydrogen content 
falls below A=0.3, i.e. when two conditions are fulfilled: 
(1) the H-burning or He-core is revealed on the surface 
(the stellar mass must have fallen below 1/2 of its initial 
value) and (2) the hydrogen content in the core, given by 
Eq. (|), is lower than X=0.3. 

In the papers referenced above about VMOs, later 
burning processes build a CO-core which almost fills 
the whole He-core (Mco^^He)- This result cannot hold 
here, since we include strong mass-loss and a possible WR 
stage which sensitively depletes the mass. The CO-core 
is assumed to fill the total remaining mass at the end of 
the lifetime. The assumed lifetimes and H-burning core 
masses, and the obtained CO-core masses, are listed in 
Tab. O. We notice that, for the higher metallicity cases, 
the final masses of the CO-cores tend to converge to the 
same small value for all the VMOs, due to a prolonged 
WR stage with very efficient mass loss. 

To derive the composition of the ejected material, we 
assume that the layers over the H-burning core have not 
altered their initial composition, i.e. we neglect any con- 
vective dredge-up of processed material and a possible re- 
cession of the convective core from initial masses larger 
than 1/2 the total mass. When the core is revealed on the 
surface, we assume that all initial CNO isotopes have been 
converted to 14 N, that the hydrogen content is given by 
Eq. (^|) at any time t and that the remaining fraction is 
4 He. 

The size of the CO-core eventually drives the SN ex- 
plosion, as we have assumed to hold for lower masses (§ |J). 
In most cases the final outcome is an iron -core collapse 
SN, whose ejecta are calculated as in § |43|. The only ex- 
ceptions are the cases M=150 and 200 M , Z=0.0004, 
whose CO-cores are marked "PC SN" in Tab. [TiJ These 
objects never lose enough mass to let the core uncovered; 



they build a 75 Mq and a 100 Mq He-core, respectively, 
which eventually explode as PC SNae. The ejecta of such 
cores are taken directly from Tab. 18 of Woosley (1986), 
but for those of 14 N which are to be scaled with metal- 
licity. In fact, Woosley's PC SNas are calculated for solar 
metallicity and their initial abundances of CNO isotopes, 
which are turned to 14 N during H-burning, are solar; for 
our Z=0.0004 case, the initial CNO content is much lower. 

The resulting ejecta for the whole set of VMOs and for 
the five metallicity cases included in our grid are listed in 
Tab. 
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8. The chemical evolution model 

Our results about stellar nucleosynthesis and ejecta de- 
scribed in previous sections were inserted in a model of 
chemical evolution, which was applied to the analysis of 
the Solar Neighbourhood. In this section we outline the 
main features and equations of the model. 

We adopt an open model with continuous infall of pri- 
mordial gas that builds the disk gradually, as suggested 
by dynamical studies (Larson 1976; Burkert et al. 1992; 
Carraro et al. 1997). Infall also provides for the better so- 
lution to the "G-dwarf problem" (Lynden-Bell 1975; Tins- 
ley 1980; Chiosi 1978, 1980; Pagel 1989; Matteucci 1991). 
Our formulation follows the one framed by Talbot & Ar- 
nett (1971, 1973, 1975) and adapted to open models by 
Chiosi (1980). This formulation, particularly suitable for 
galactic discs, is also currently used by other authors (e.g. 
Matteucci & Francois 1989; Timmes et al. 1995). 

The Galactic disc is divided into concentric cylindri- 
cal shells, 2 kpc wide each, which evolve independently, 
neglecting any possible radial flows of gas or stars. Each 
ring consists of a homogeneous mixture of gas and stars, 
so that for each ring the only independent variable is time 
t (one-zone formulation, Talbot & Arnett 1971). Galactic 
discs are comfortably described in terms of surface mass 
density o~(r,t), which depends both on the galactocentric 
radius r of the ring and also on time t, since in each ring 
the surface density is growing in time due to gradual infall 
of gas. Since in this paper the model is applied only to the 
Solar Neighbourhood (r = r Q ), for the sake of simplicity 
we drop the dependence on r in the following equations. 
If we indicate with a g (t) the surface gas density, the gas 
fraction in the ring at any time t is the ratio: 



while the surface star density is: 
a s (t) = a(t) - a g (t) 

In closed models the total surface density a is constant and 
the other quantities can be normalized with respect to it; 
in open models it is suitable to normalize with respect to 
the total surface density at the present age of the Galaxy 
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12. Total ejecta in Mq for VMOs (150 -f- 1000 M©) of different metallicities 



Z=0.0004 



M 


L H 




4 Ho 


V2 C 


is c 


i4 N 




16 O 


iV o 


1H 












*Nc 


' z ' z Nc 


M Mg 


*Si 


sa s 


4U Ca 


bb Fo 


M r 


150 


59.51 


6.60E-05 


20.95 


1.31 


8.88E-05 


2.67E-03 


7.32E-06 


38.02 


6.30E-06 


3.63E-05 












1.80 


3.08E-04 


1.60 


16.00 


8.90 


1.30 


0.612 





200 


79.82 


8.80E-05 


28.90 


0.840 


1.18E-04 


1.46E-02 


9.76E-06 


34.02 


8.40E-06 


4.84E-05 












1.50 


4.10E-04 


1.50 


23.00 


15.00 


2.40 


13.00 





300 


152.7 


1.32E-04 


129.1 


1.48E-02 


1.78E-04 


4.48E-02 


1.46E-05 


3.17E-02 


1.26E-05 


7.27E-05 












8.61E-03 


1.16E-03 


2.66E-03 


3.68E-03 


2.23E-03 


3.38E-04 


6.60E-03 


18.00 


500 


292.2 


2.20E-04 


186.6 


2.47E-02 


2.96E-04 


7.75E-02 


2.44E-05 


5.28E-02 


2.10E-05 


1.21E-04 












1.46E-02 


1.97E-03 


4.52E-03 


6.26E-03 


3.79E-03 


5.74E-04 


1.12E-02 


21.00 


1000 


596.2 


4.40E-04 


382.4 


4.94E-02 


5.92E-04 


0.162 


4.88E-05 


0.106 


4.20E-05 


2.42E-04 












2.99E-02 


4.02E-03 


9.24E-03 


1.28E-02 


7.75E-03 


1.17E-03 


2.29E-02 


21.00 


Z=0.004 
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4 Hc 


C 


13 c 
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17 O 


1S Q 
















^Mg 




32 g 


4U Ca 


! *Fc 


M r 
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56.67 


6.60E-04 


57.84 


7.41E-02 


8.88E-04 


0.153 


7.32E-05 


0.159 


6.30E-05 


3.63E-04 












3.51E-02 


4.72E-03 


1.09E-02 


1.50E-02 


9.11E-03 


1.38E-03 


2.69E-02 


35.00 


200 


75.57 


8.80E-04 


88.74 


9.89E-02 


1.18E-03 


0.227 


9.76E-05 


0.211 


8.40E-05 


4.84E-04 












5.04E-02 


6.77E-03 


1.56E-02 


2.15E-02 


1.31E-02 


1.98E-03 


3.86E-02 


35.00 


300 


159.2 


1.32E-03 


131.6 


0.442 


1.78E-03 


0.480 


2.73E-04 


4.34 


1.26E-04 


7.27E-04 












0.692 


1.20E-02 


0.330 


0.126 


3.01E-02 


4.10E-03 


7.29E-02 


2.68 


500 


295.9 


2.20E-03 


192.0 


0.551 


2.96E-03 


0.810 


3.18E-04 


3.82 


2.10E-04 


1.21E-03 












0.658 


2.01E-02 


0.226 


9.95E-02 


4.01E-02 


5.92E-03 


0.115 


5.67 


1000 


591.9 


4.40E-03 


394.0 


0.798 


5.92E-03 


1.65 


5.62E-04 


4.35 


4.20E-04 


2.42E-03 












0.811 


4.02E-02 


0.273 


0.165 


7.97E-02 


1.19E-02 


0.232 


5.67 


Z=0.008 
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^Mg 
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4U Ca 
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M r 
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55.54 


1.32E-03 


80.25 


0.460 


1.78E-03 


0.423 


1.75E-04 


2.86 


1.26E-04 


7.27E-04 












0.636 


1.12E-02 


9.16E-02 


4.42E-02 


2.22E-02 


3.32E-03 


6.41E-02 


9.57 


200 


80.98 


1.76E-03 


104.4 


0.509 


2.37E-03 


0.592 


2.24E-04 


2.97 


1.68E-04 


9.69E-04 












0.667 


1.53E-02 


0.101 


5.73E-02 


3.02E-02 


4.52E-03 


8.75E-02 


9.57 


300 


159.1 


2.64E-03 


131.9 


0.572 


3.35E-03 


0.968 


4.31E-04 


3.81 


2.52E-04 


1.45E-03 












0.433 


2.41E-02 


0.158 


0.383 


0.197 


2.28E-02 


0.279 


2.15 


500 


292.9 


4.40E-03 


195.0 


0.778 


5.92E-03 


1.63 


6.13E-04 


4.92 


4.20E-04 


2.42E-03 












0.841 


4.04E-02 


0.361 


0.227 


9.48E-02 


1.35E-02 


0.244 


2.96 


1000 


582.7 


8.80E-03 


401.1 


1.27 


1.18E-02 


3.31 


1.10E-03 


5.97 


8.40E-04 


4.84E-03 












1.14 


8.14E-02 


0.456 


0.358 


0.174 


2.55E-02 


0.478 


2.96 
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^Nc 


" 4 Mg 


*Si 
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M r 
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52.39 


3.30E-03 


88.62 


0.614 


4.44E-03 


1.17 


5.04E-04 


3.50 


3.15E-04 


1.82E-03 












0.429 


2.95E-02 


0.134 


0.417 


0.206 


2.46E-02 


0.334 


2.11 


200 


96.71 


4.40E-03 


94.51 


0.672 


5.92E-03 


1.61 


5.73E-04 


2.89 


4.20E-04 


2.42E-03 












0.379 


4.01E-02 


0.116 


0.406 


0.223 


2.51E-02 


0.317 


2.09 


300 


155.3 


6.60E-03 


133.4 


0.919 


8.88E-03 


2.45 


8.17E-04 


3.42 


6.30E-04 


3.63E-03 












0.531 


6.06E-02 


0.613 


0.471 


0.263 


3.11E-02 


0.434 


2.09 


500 


281.4 


1.10E-02 


202.9 


1.44 


1.48E-02 


4.12 


1.33E-03 


4.90 


1.05E-03 


6.05E-03 












0.799 


0.101 


0.266 


0.668 


0.361 


4.59E-02 


0.784 


2.09 


1000 


553.8 


2.20E-02 


420.2 


2.67 


2.96E-02 


8.32 


2.55E-03 


7.54 


2.10E-03 


1.21E-02 












1.56 


0.204 


0.502 


0.994 


0.559 


7.59E-02 


1.37 


2.09 


Z=0.05 
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a He 


4 Ho 
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i« 
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^Nc 


M M S 


28g ; 


32 g 


4U Ca 


bB Fo 


M r 


150 


44.58 


8.25E-03 


94.29 


1.14 


1.11E-02 


3.02 


1.01E-03 


3.00 


7.87E-04 


4.54E-03 












0.813 


7.51E-02 


0.221 


0.367 


0.194 


2.79E-02 


0.494 


1.74 


200 


87.44 


1.10E-02 


98.84 


1.45 


1.48E-02 


4.07 


1.31E-03 


3.66 


1.05E-03 


6.06E-03 












1.00 


0.101 


0.280 


0.449 


0.243 


3.54E-02 


0.641 


1.74 


300 


140.9 


1.65E-02 


140.2 


2.07 


2.22E-02 


6.17 


1.92E-03 


4.98 


1.58E-03 


9.08E-03 












1.39 


0.152 


0.398 


0.612 


0.342 


5.03E-02 


0.933 


1.74 


500 


246.6 


2.75E-02 


224.1 


3.30 


3.70E-02 


10.37 


3.14E-03 


7.62 


2.62E-03 


1.51E-02 












2.15 


0.255 


0.634 


0.938 


0.540 


8.03E-02 


1.52 


1.74 


1000 


478.2 


5.50E-02 


466.6 


6.39 


7.40E-02 


20.88 


6.19E-03 


14.23 


5.25E-03 


3.03E-02 












4.06 


0.511 


1.22 


1.75 


1.04 


0.155 


2.98 


1.74 
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tc i.e. at the final age of the model. So we introduce the The contribution of the infalling gas to the evolution of 
normalized surface gas density as: the (normalized) gas fraction G(t) is: 



G(t) 



dt 



G(t) 



inf 



*(*G) 



(11) 



In each ring the gas is assumed to be chemically homoge- 

ncous, and the normalized gas density for each chemical and < for each sin § lc chcmical s P ecies 
species i is: 



G i (t)=X i (t)G(t) 

where X.- L is the fractionary mass abundance of species i; 
J2i Xi=l by definition. The chemical evolution of the ISM 
is the evolution of the set of the Gj's, described by: 



,inf 



inf 



which gives the third term on the right side of Eq. (7). 



(12) 



+ *(t - t m ) Rm(t - tmMM) dM + 

JM, 

(7) 



8.2. The Star Formation Rate 

We adopt the formulation of the SFR for the Galactic disk 
given by Talbot & Arnett (1975): 



+ 



M, 



d , . 



<r(r,t) a g {r, t) 
a(f,t) 



,{r,t) 



(13) 



inf 



where *(f) is the SFR, $(M) is the IMF, R m (t) is the 
mass fraction of a star of mass M ejected into the ISM 
in the form of element i, Mi and M u are the lower and 
upper limit for stellar mass respectively and tm is the 
lifetime of a star of mass M. The first term on the right 
side represents the depletion of species i from the ISM due 
to star formation; the second term represents the amount 
of species i put in the ISM by stellar ejecta; the third 
term is the contribution of the infalling gas. We describe 
the various ingredients in the following subsections. 

8.1. The inf all term 

In open models the surface mass density o~(r,t) increases 
by slowly accreting gas at a rate a in f (r, t) , until it reaches 
the observed present values. An infall rate exponentially 
decreasing in time with a timescale r: 



2f (r,t) = A{r)e- 



(8) 



Here, k is the exponent of the Schmidt (1959) law for star 
formation, ^ oc p K ; plausible values for k range from 1 to 
2 (proportional to gas density or proportional to cloud- 
cloud collision events respectively, Larson 1991); a(r,t) is 
the surface mass density at a given galactocentric distance 
f, adopted as a normalization factor; v is a parameter for 
the star formation efficiency, related to the choice of f . 

Talbot & Arnett's SFR is based on Schmidt's law, but 
it also takes into account that the cooling due to gas accre- 
tion onto the equatorial plane is balanced by the heating 
due to the feed-back of massive stars. The SFR turns out 
to be related to the dynamical timescale, shorter where the 
mass density is larger, while the timescale of gas accretion 
onto the equatorial plane is longer than the purely dy- 
namical one due to the feed-back from star formation. As 
a caveat, we remind that the formula above holds in the 
plane-parallel approximation, which may not be a good 
one in the case of open models. 

In terms of the formalism introduced in previous sec- 
tions, this adopted SFR 



G K (r,t) (14) 



well reproduces the results of dynamical models (Larson 
1976; Burkcrt et al. 1992; Carraro et al. 1997), with the *M) = v 
exception of radial flows. A(r) is obtained by integrating 

upon time and by imposing that at the age t G the observed Limiting to the solar ring and taking f = r& we get (drop _ 



~a(r,ty 


2(k-1) 


'a(r, t G )' 


K-l 


[a{f, t)_ 




[ <r(r,t) _ 





present surface mass density a(r,tc) is matched: 
A(r) (l-e - ^) T = <r(r,t G ) 
In the case of the solar ring: 
A Q (l - e~"^) t = cr(r Q ,i G ) 



ping the dependence on r): 



(9) 



*(*) 



<r(*c) 
a(t) 



G K {t) 



(15) 



where v is to be fixed so as to reproduce the features of 

(10) 

the Solar Neighbourhood. 
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Table 13. The Qij matrix 





(1) 


4 ff e 

(2) 


VI c 

(3) 


(4) 


14 N 

(S) 


lb 

(6) 


(7) (8) 


24 M „ 

(9) 


2 8 Si 
(10) 


32 s 
(11) 


4U Ca 

(12) 


bb Fc 
(13) 


(1) 


H 


1 — 94 
























(2) 
(3) 


*He 
12 C 


XC™C 


1 - q c 

XC™C 


1 - "!C13s 




















(4) 


13 C 


*C13™C 


XC13™C 


9C13s - INs 


1 - INs 


















(S) 


14 „ 


XjV™C 


XJV™C 


INs - 1C 


INs - 1C 


1 - 1C 


INs - 1C 














(6) 


16 Q 


xo™C 


xo™C 








1 - l»i 














(7) 

(8) 


nr 

20]V e 


M«"C 


XNe^'C 


VJQ 


WQ 






1 - d 

1 - d 












(9) 


24 Mg 


XMg w C 


XMg™C 












1 - d 










(10) 


28 Si 


XSi™C 


XSi™C 














1 - d 








(11) 


32 S 


XS«"C 


XS«"C 
















1 - d 






(12) 


«C» 


XCa™C 


XCa™C 


















1 - d 




(13) 


56 Fe 


XF„»r7 






















1 - d 



8.3. The Initial Mass Function 
We adopt a Salpeter-like IMF: 

$(M)dM = CM-^ dM 



(16) 



where /Lt=1.35 in Salpeter's (1955) law. 

Once the slope \x and the limiting masses Mi and M u 
are chosen, the coefficient C is fixed by normalizing the 
IMF over the whole mass interval: 



M„ 



$(M) dM = 1 



(17) 



Mi 



Since the bulk of chemical enrichment is due to stars with 
M >1 M Q , it is meaningful to fix the fraction £ of the 
total stellar mass distributed in stars above 1 M , which 
is equivalent to fixing Mi (see below). Then we get the 
normalization condition: 



<$>{M)dM = C / M'^dM 



(18) 



where M x =l M G and M u is chosen to be 100 M Q (no 
VMOs are supposed to be present). The normalization 
condition fixes Mi once £ is given: 



/■Mi 

/ $(M)dM + C 



(19) 



The slope /z of the IMF may not be constant over the whole 
range of stellar masses (Miller & Scalo 1980; Scalo 1986). 
The model includes also the possibility of a variable \x 
over different mass ranges; in this case, beside the normal- 
ization condition, we impose that the IMF is continuous 
where the different mass ranges connect and thus deter- 
mine the normalization coefficients Ci, C n for the dif- 
ferent ranges. See, for instance, Chiosi & Matteucci (1982). 



8.4- The contribution of stellar ejecta 

We calculate the fraction of a star of initial mass M that 
is ejected back in form of chemical species i 



Rmi — 



EiM 

1m~ 



(20) 



according to our results on stellar ejecta (discussed in pre- 
vious sections) and by means of the "Qij matrix" formal- 
ism, first introduced by Talbot & Arnett (1973) and later 
adopted by many authors (Chiosi & Matteucci 1982; Mat- 
teucci & Frangois 1989; Ferrini et al. 1992). Each matrix 
element Qij is defined as the mass fraction of a star orig- 
inally in form of species j which has been processed and 
ejected as species i. Rmi is the sum over all "fuels" j: 



Rmi 



N 

Eg*. 

3=1 



(M) Xj 



Each matrix clement is defined as: 



Qij 



M, 



XjM 



(21) 



(22) 



where Xj is the initial mass abundance of species j and 



Mi 



is the amount of species i synthesized starting from 



j and eventually expelled. The diagonal elements Qu rep- 
resent the unprocessed fraction of species i that is eventu- 
ally re-ejected. By summing upon all the sources of species 
i, one must get the total ejecta of i: 



EQ 



i 3 XjM = EiM 



(23) 



Each stellar mass M corresponds to a different Qij matrix, 
which generally depends also on the metallicity of the star. 

The Qij matrix of Talbot & Arnett originally treated 
!H, 2 H, 3 Hc, 4 He, 12 C- 16 0, 14 N, neutron rich isotopes {nr) 
and the bulk of heavy elements (h). Later, the dimensions 
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of the matrix were extended to include a higher number 
of chemical species (Ferrini et al. 1992). The formalism we 
use here is analogous to that in Ferrini et al. (1992), al- 
though the definitions of some matrix elements have been 
revised In our model we follow the evolution of 13 ele- 
ments: 1 H, 4 He, 12 C, 13 C, 14 N, 16 0, neutron rich isotopes 
(nr), 20 Ne, 24 Mg, 28 Si, 32 S, 40 Ca, 56 Fe. 

The set of non-zero Qij elements is displayed in 
Tab. |l3]. Here we list the various quantities entering the 
matrix elements: 

d mass fraction eventually locked in 

the remnant 
54 mass fraction involved in H-burning 

qc mass fraction involved in He- 

burning 

w c — Qc ~ d mass fraction that has been pro- 
cessed by He-burning and has been 
ejected 

qNs mass fraction where 12 C, 13 C and 

16 O have been turned to secondary 
14 N by the CNO cycle 
qci3s mass fraction where 12 C is turned to 

secondary 13 C by the CNO cycle 
Xi fractionary mass abundance within 

wc of newly synthesized species i 

The meaning and the calculation of the non-zero elements 
of the Qij matrix are discussed in App. C. 

8.5. Stellar lifetimes 

To follow the temporal behaviour of different chemical 
species and isotopes, we drop the "instantaneous recycling 
approximation" by taking into account the role of finite 
lifetimes for stars of different masses. Our stellar models 
of different metallicities allow us to consider the effects of 
a different metal content not only on the ejecta but also 
on the lifetimes tm = tm{Z). The dependence of lifetimes 
on metallicity can be up to a factor of 2 for stars around 
1M . 

The lifetimes we adopt are calculated as the sum 
t = tu + tue of the H-burning and He-burning timescales 
of the stellar tracks of the Padua library, and are listed in 
Tab. [li]. The trend of tm with metallicity is not univocal: 
up to solar metallicity, tm increases with Z due to the 
effect of increasing opacity, while for super-solar metal- 
licity (Z =0.05) tm decreases because the helium con- 
tent Y increases, due to the assumed ratio AY/ AZ=2.5. 
When Y increases sensitively, (1) the hydrogen content 
X = 1 — Y — Z, i.e. the available "fuel" , decreases and (2) 
the average molecular weight fj, increases, which leads to a 
higher luminosity (L ~ M 7 ' 4 )- For super-solar metallicities 
the effects of Y overcome those of opacity, giving shorter 
lifetimes with increasing Z . 

In the chemical model, for each mass M and metallic- 
ity Z we calculate the corresponding lifetime tm(Z) by in- 
terpolating within the logarithmic relation log(M) — log(i) 



for the tabulated metallicities, and then by interpolating 
with respect to Z. 

Around ~1 M the influence of metallicity reduces 
lifetimes of a factor of 2 for low metallicities with respect 
to the solar case: low and intermediate mass stars born in 
early galactic stages eject their nucleosynthetic products 
in shorter times. 

We assume that each star expels its ejecta all at once 
at the end of its lifetime, and that the ejected material is 
immediately mixed in the ISM, which remains always ho- 
mogeneous. This "instantaneous mixing approximation" 
is suitable to reproduce the average trends of the age- 
metallicity relation, of abundance ratios and so on, while 
it can't model the observed scatter of the data around 
the average trend. Only few models relaxing the instan- 
taneous mixing approximation can be found in literature 
(Malinie et al. 1991, 1993; Pilyugin & Edmunds 1996; van 
den Hoek & de Jong 1997). 

8.6. The numerical solution 

For the set of equations 7 we adopt the numerical solution 
suggested by Talbot & Arnett (1971), with the addition 
of the infall term (Chiosi 1980). Each equation is treated 
as a linear differential equation, whose analytical solution 
can be written as: 

Gi(t 2 ) = Giih) e -x(^) + f 2 Witye-xW dt (24) 
where: 



X (h,t 2 )= / ri(t)dt, 



_ B(t) _ $(*) 
V[ > ~ G(t) ~ G(t) 



and 

Wi(t) = 



M„ 



* (t - T M )HM)R Mi {t - T M ) dM 



M, 



dt 



Gi(t) 



(25) 



We perform the integration over a timestep At — t n+1 —t n 
by approximating 77 and Wi over the whole timestep with 
a constant value (the one estimated at time t n+ 2 = t n + 
lAi): 



G l (t n+1 ) = G t (t n )e^ At + — [1 - e-" At l 

r\ J 



(26) 



Since rj and Wj are not strictly constant over At, an itera- 
tion is required to provide with convergency. As suggested 
by Talbot & Arnett (1971), we choose to iterate only with 
respect to 77, because Wi is an integrated quantity involv- 
ing all past value of G(t) and is not so sensitive to the 
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Table 14. Stellar lifetimes in years for stars from 0.6 to 120 Mq, for different metallicities 



M 


Z=0.0004 


Z=0.004 


Z=0.008 


Z=0.02 


Z=0.05 


0.6 


4.28E+10 


5.35E+10 


6.47E+10 


7.92E+10 


7.18E+10 


0.7 


2.37E+10 


2.95E+10 


3.54E+10 


4.45E+10 


4.00E+10 


0.8 


1.41E+10 


1.73E+10 


2.09E+10 


2.61E+10 


2.33E+10 


0.9 


8.97E+09 


1.09E+10 


1.30E+10 


1.59E+10 


1.42E+10 


1.0 


6.03E+09 


7.13E+09 


8.46E+09 


1.03E+10 


8.88E+09 


1.1 


4.23E+09 


4.93E+09 


5.72E+09 


6.89E+09 


5.95E+09 


1.2 


3.08E+09 


3.52E+09 


4.12E+09 


4.73E+09 


4.39E+09 


1.3 


2.34E+09 


2.64E+09 


2.92E+09 


3.59E+09 


3.37E+09 


1.4 


1.92E+09 


2.39E+09 


2.36E+09 


2.87E+09 


3.10E+09 


1.5 


1.66E+09 


1.95E+09 


2.18E+09 


2.64E+09 


2.51E+09 


1.6 


1.39E+09 


1.63E+09 


1.82E+09 


2.18E+09 


2.06E+09 


1.7 


1.18E+09 


1.28E+09 


1.58E+09 


1.84E+09 


1.76E+09 


1.8 


1.11E+09 


1.25E+09 


1.41E+09 


1.59E+09 


1.51E+09 


1.9 


9.66E+08 


1.23E+09 


1.25E+09 


1.38E+09 


1.34E+09 


2.0 


8.33E+08 


1.08E+09 


1.23E+09 


1.21E+09 


1.24E+09 


2.5 


4.64E+08 


5.98E+08 


6.86E+08 


7.64E+08 


6.58E+08 


3 


3.03E+08 


3.67E+08 


4.12E+08 


4.56E+08 


3.81E+08 


4 


1.61E+08 


1.82E+08 


1.93E+08 


2.03E+08 


1.64E+08 


5 


1.01E+08 


1.11E+08 


1.15E+08 


1.15E+08 


8.91E+07 


6 


7.15E+07 


7.62E+07 


7.71E+07 


7.45E+07 


5.67E+07 


7 


5.33E+07 


5.61E+07 


5.59E+07 


5.31E+07 


3.97E+07 


9 


3.42E+07 


3.51E+07 


3.44E+07 


3.17E+07 


2.33E+07 


12 


2.13E+07 


2.14E+07 


2.10E+07 


1.89E+07 


1.39E+07 


15 


1.54E+07 


1.52E+07 


1.49E+07 


1.33E+07 


9.95E+06 


20 


1.06E+07 


1.05E+07 


1.01E+07 


9.15E+06 


6.99E+06 


30 


6.90E+06 


6.85E+06 


6.65E+06 


6.13E+06 


5.15E+06 


40 


5.45E+06 


5.44E+06 


5.30E+06 


5.12E+06 


4.34E+06 


60 


4.20E+06 


4.19E+06 


4.15E+06 


4.12E+06 


3.62E+06 


100 


3.32E+06 


3.38E+06 


3.44E+06 


3.39E+06 


3.11E+06 


120 


3.11E+06 


3.23E+06 


3.32E+06 


3.23E+06 


3.11E+06 



exact value of G(t n+1 ), and also because the resolving ex- 
pression depends exponentially on 77, while only linearly 
on Wi. If Gf ] is the estimate of G l {t n+1 ) at the fcth iter- 
ation, we proceed through successive corrections: 



(fe+i) = G (fc) 



G 



where 



1 



sr(fc+i) 



(fc+l) _ 



A - G 



G t (t n )-^ 



1 d In T] 

lil= 2d^G X 



w 

V 



(27) 



(28) 



(29) 



In this way, the infall term in the Wi 's can be assumed to 
be constant within At, and the formulation by Talbot & 
Arnett (1971), suited to close models, can be applied also 
to open models. 

As for the initial conditions, we assume c(0) to be 
very small, although non-zero to avoid mathematical in- 
finities for t=0. At the beginning the disk is formed by 
gas only, while star formation is not active yet: cr g (0) = 
er(0). With the adopted normalization, this translates in 
G(0) = G,(0) = G s (0)~0. 

The Wi& are calculated by integrating with respect 
to time, rather than with respect to mass; Eq. (25) is re- 
written as: 



Wi(t) 



t — TM u 



<S>(M)R Mi (t') 



dM 



dt'- 



M(t-i') 



,- v At 



V 



rjAte-^ \ Gi(t n ) 



(30) 



and 77, Wi are estimated at t n+ 2 . G\ is updated until 
the correction 8\ k+1 ^ is smaller than a chosen limit 8 max - 
The timestep At of the model is chosen to be the min- 
imum value between: (1) a timestep Ati which guarantees 
that the relative variation of the Gi(?;)'s is lower than a 
fixed e; (2) a timestep Ai 2 which guarantees that the sur- 
face mass density er(i) increases no more than 5%; (3) a 
timestep At3 which is twice the previous timestep of the 
model, in order to speed up the calculation when possible. 



dt 



Giit) 



(31) 



where M(t) is the mass of a star of lifetime r. We need 
to integrate with respect to time when introducing the 
dependence of metallicity, since all quantities then depend 
on Z(t) as well as on M: the integral on the right hand 
is to be calculated on the path of the (M, t) plane fixed 
by the relation between mass and lifetime, or better by 
the M(r,Z(t — r)) vs. r relation, where Z{t) is built by 
the on-going model its elf. Also the "restitution fractions" 
Rm% depend on Z{t) (§ |S.4|) ; therefore, all the stored values 
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of G l (t n ), Xi{t n ) and Z{t n ) of all the timesteps t n of the 
model enter the evaluation of ^(f), B,Mi(t') and dM/dTM 
for each t' in the integral. 

Operatively, the integration is turned to a summation 
over a series of time intervals [#-!),#)], with = 
e f( femt ") = t (the current age of the model). Within each 
time interval, the integrand is assumed to be constant and 
is estimated at the middle point (i^ 5-1 ) +t^)/2. The time 
intervals [#-!),#)] arc fixed as follows. The whole mass 
range [M U M U ] is divided in 200 intervals [JlfW,!^ 1 )], 
where the MW are equidistant in logarithmic scale. At 
any age t of the model, we calculate the corresponding 
birthtimes = t — t m m . The 's defining the grid for 
the integration are given by the t^ p \ implemented with 
the previous timesteps t n of the model. 

8.7. Inserting Type la supernova; 

Up to now we have discussed the ejecta of single stars of 
different mass ranges, but in a chemical model we also 
need to include Type la supernovae (SNae la), that origi- 
nate in close binary systems and contribute an important 
fraction of heavy elements, especially iron. We need to 
introduce SNae la especially to explain the observed evo- 
lution of abundance ratios of a-elemcnts with respect to 
iron (Sncdcn ct al. 1979; Matteucci & Greggio 1986). 

Prescriptions for the rate of SNae la and for the com- 
position of their ejecta are needed. We adopt the rate sug- 
gested by Greggio & Renzini (1983), which assumes the 
scenario of Whelan & Iben (1973): SNae la, are due to 
the explosion of a CO white dwarf that reaches the Chan- 
drasekhar limit by accreting material from a giant com- 
panion filling its Roche lobe. An upper limit to the mass 
of the primary is fixed by the requirement that it builds a 
degenerate CO-core before filling its Roche lobe, so that 
it doesn't explode as a SN II. For models with convective 
overshooting, this translates in 

Mi < 6Mq 

where Mi is the mass of the primary (Mi > M 2 ); there- 
fore, the upper limit for the total mass Mb = Mi + M 2 
of a binary system able to produce a SN la is: 



M 



B,u 



\2M P 



© 



Since the primary needs to accrete enough mass to 
reach the Chandrasekhar mass, a minimum total mass 
for the system is also introduced, generally Mb,i ~3 M Q . 
The distribution function of the fractionary mass of the 
secondary [i — M 2 /Mb, ^ < 1/2 is (Greggio & Renzini 
1983): 

F( M ) - 2 V 

which is normalized between and 0.5 and favours sys- 
tems whose components have similar masses (Mi ~ M 2 , 
Tutukov & Yungelson 1980). 



As the star formation rate for single stars of mass M 
is \E'(t) < I>(M), we assume that the SFR for binary system 
precursors of SNae la is A &(t)&(M). This means that a 
fraction A of stars between Mb,i and Mb,« is assumed 
to form binaries (with such characteristics that a SN la 
is eventually produced), rather than single stars. A is a 
parameter to be fixed so to match the observed rate of 
SNae la. 

In this scenario, the typical timescale for the explosion 
of a SN la is fixed by the lifetime tm 2 01 the secondary. 
The explosion rate (by number) of SNae la is expressed as 
(Greggio & Renzini 1983): 



R S Ni(t) 



/ ^r 1 / m*(t-T Ma )diM 



(1Mb (32) 



where /i m is the minimum mass fraction contributing to 
the rate of SNae la at time t: 



ji m = max 



M 2 (t) M B -0.5M B . 



Mf 



For Type la SNae we adopt the ejecta Esni of the W7 
model of Nomoto et al. (1984), in the updated version by 
Thielemann et al. (1993). Due to the homogeneity of SNae 
la, we can safely assume that the same set of ejecta holds 
for all of these objects. 

To include the contribution of SNae la, we need to 
change the formulation of the Wi's. Following Matteucci 
& Greggio (1986) with some slight changes, in the range 
Mb.i — Mb, u we distinguish the contribution of single stars 
- a fraction (1 — A) of the whole — from that of bina- 
ries originating SNae la — a fraction A of the whole. In 
the case of the binaries, we assume that their ejecta are 
released in two steps: after a time tm 1 the primary expels 
its products behaving just like a single star (i.e. according 
to the Qij matrix suited to its mass and metallicity) , while 
after a time tm 2 the secondary pours mass on the compan- 
ion originating the SN explosion. Therefore, introducing 
the notation: 

^R^M^t) = *(i - T M MM)R Mi {t - t m ) 
we can write: 



Wi(t) = / ^^>R l (M,t)dM - 

I Mi 



(1-A) / *$i?,(M,i)(iM- 
'm b . 



f 

/ ■$<f>R i (M,t)dM + 

J M BjU 



r M B,u <j)f M R ) 
- A I -A-^ [h (t, M B ) + h (t, M B )] dM B + 

M B j M B 
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(33) 



in} 



indicating with: 



h(t, M B ) = / F(ti)V{t - tmJRmA* ~ T Ml )M!dfi 



- T M2 ) EsNIidfl, 



Ah = (1 - fj)M B 
h(t,M B ) = 
M 2 = (iM B 



In Eq. (|33|), the first three terms on the right hand side 
represent the contribution of single stars, the fourth term 
represents that of binaries becoming SNae la and the fifth 
term is the contribution of infall. Since the ejecta Eg^n 
of SNae la are assumed to be independent of M B or /i, the 
term describing binaries can also be written as: 



A / \/> h (t, M B ) (1Mb + RsNi(t)E SN ii 

JMba M B 



This formulation better shows how the primary star is as- 
sumed to evolve like a single star, unaffected by its com- 
panion as far as nucleosynthesis is concerned and releasing 
all its ejecta away from the system, while the secondary 
pours all its ejecta to produce the SN la, without any 
direct contribution to the enrichment of the ISM. 

The term expressing the ejecta of the primary stars can 
also be integrated with respect to the mass of the primary 
Mi, if we write it as: 

rM ltmax 

A *(t - t Mi ) R Mll (t - r{Ah) JF(Mi) dM x 



where: 



Mi 



^(Mi) 



M 



B.l 



AI 



l.max 



M 



B.u 



/(//)<!> du 



Mx 



f{v) = 2A{l-v)\ ue [0.5,1] 



max < 0.5, 



Mi 



M 



B.u 



min < 1 



Mi 

w 



B.l 



Integrating the W% 's with respect to time and introducing 
the notation: 



$RidM(t,t') = 
TR l dM 1 (t,t') = 



, , dM\ 



M(i-t') 



F(Mi)R Mli tf) (- 



dM x 

~dV~ 



Ah 



when we insert SNae la Eq. becomes 



Wi(t) = 



t — TMp 







$RidM(t,t')dt' 



+ (1-A) V(t')®RidM(t,t')dt' + 

Jt-^B.l 

pt—T M u 

+ y(t')<f>RidM(t,t')dt' + 

'~ TM B,« 

rt— TMi mam 

+ A ^{t')J 7 R l (t')dM 1 (t,t')dt / + 

J t-T Ml min 

+RsNlEsNIi + 
d 



dt 



Gi(t) 



(34) 



The integration is performed as described in § 8.6. Again, 
here we need to integrate with respect to time rather than 
with respect to mass if we want to include the implicit 
dependence on the metallicity Z(t). 

9. Chemical evolution in the Solar Neighbourhood 

As a first application of our model, we analyze the chemi- 
cal evolution in the Solar Neighbourhood. Standard obser- 
vational counterparts for chemical models, such as (1) the 
current gas fraction, (2) the rate of Type I and Type II 
SNae (3) the age-metallicity relation, (4) the past and 
current estimated SFR, (5) the distribution of long-lived 
stars in metallicity (G-dwarf problem), are used to cali- 
brate the free parameters of the model. The abundance 
ratios of different elements observed in the atmospheres 
of nearby stars are used as a test for our nucleosynthesis 
prescriptions, since model predictions for most abundance 
ratios depend mainly on the adopted yields and IMF, and 
only slightly on other model parameters. 
Here we list the parameters of the model: 



V 


star formation efficiency 


K 


exponent of the star formation law 


M 


exponent of the IMF 


M u 


higher limit for the IMF 


C 


mass fraction of the IMF in stars with 




M >1 M Q 


T 


infall timescale 


to 


age of the system 


A 


"amplitude" factor for SNae la 




galactoccntric radius of the Sun 


o"O* ) 


surface mass density in the Solar 




Neighbourhood 



Mi(t-t') 



Some of these parameters are directly fixed by observa- 
tional determinations. We adopt r ~ 8 kpc, as indicated 
by recent estimates of the galactocentric distance of the 
Sun (Reid 1993; Paczynski & Stanek 1997). From a col- 
lection of independent measurements of the local surface 



Portinari et al. : Chemical evolution with metallicity dependent yields 



23 



mass density of the disc Sackctt (1997) derives a fiducial 
estimate of 53±13 M©/pc 2 ; in this paper we therefore as- 
sume er(r Q )=50 M©/pc 2 . 

For the local IMF, we assume a Salpeter (1955) slope 
^=1.35 for M <2 M Q and, following Scalo (1986), /x=1.7 
for M >2 Mq. Actually, the assumed slope below 1 M Q 
is relatively unimportant for our models, since these stars 
give little contribution to the chemical enrichment and 
only act as a sink of matter from the gaseous phase. What 
actually matters in the model is the fraction of mass that 
each stellar generation stores in such low-mass, long lived 
stars, which is formulated in term of a free parameter as 

(K). 

We further assume the infalling material to be pro- 
togalactic gas with primordial chemical composition: 76% 
*H and 24% 4 He (Walker et al. 1991); the abundance of all 
heavier elements is assumed to be zero in the infalling gas 
(in practice, we assume a negligible abundance of 1CP 11 
to avoid mathematical infinities at the beginning of the 
model). 

On the other hand, v, n, r, £ and A are treated as 
free parameters to be calibrated in order to match the 
observational constraints. In the following we discuss the 
space of model parameter with respect to the available 
empirical constraints. 




Fig. 10. Current normalized infall rate vs. infall timescale, 
as predicted by model equations (see text). The shaded area 
traces the observational estimates. 



9.1. The current infall rate 

In principle, an observational estimate of the current infall 
rate on the galactic disc could fix the infall timescale r, 
since combining Eqs. (^J) and Eq. (|l^) we get 



(to) 



(.*-!) 



independently of other model parameters. 

Unluckily, current data cannot provide tight con- 
straints in practice. Observational evidence of infall comes 
from High Velocity Clouds (Oort 1970) and Very High 
Velocity Clouds, infalling on the disk at a rate of ~0.2 
and 1.0 M Q /yr, respectively (Tosi 1988 and references 
therein). While there is general agreement on the pri- 
mordial nature of VHVCs, the metal content displayed 
by some HVCs suggests a Galactic origin for at least 
a part of them (Schwarz et al. 1995). So, we can take 
the numbers above as the uncertainty range for the cur- 
rent global infall rate on the galactic disc. Assuming that 
such infall is uniform over the disc and that the radius 
of the disc is ^15 kpc, we get a local infall rate of 0.3- 
1.5 Af /pc 2 /Gyr, or & inf /a(r Q , t G ) ~ 6 HT 3 - 3 lO" 2 in 
the Solar Neighbourhood, taking er(r Q )=50 M Q /pc 2 . As 
displayed in Fig. |Io|, any r > 4 (or even less since our 
estimate is rather crude) is basically allowed within the 
uncertainties, and r is not properly constrained. 

Better hints about the infall timescale come from the 
metallicity distribution of local long-lived stars: recent 
data indicate infall timescales longer that ^5 Gyrs (see 



§ 9.6). Therefore, we start our investigation of the param- 
eter space by selecting the range r = 5 -f- 9 Gyr. 



(35) 9.2. The current gas fraction 



The current gas fraction predicted by our models is to 
be compared to observational estimates. From a compila- 
tion of data on the HI and H2 distribution, Rana (1991) 
quotes as the current surface gas density in the Solar 
Neighbourhood ^(r©) ~5.7-7 M Q /pc 2 , which implies a 
gas fraction ^0.1-0.15 of the total surface mass density 
(~5OM /pc 2 ). 

The gas fraction G(to) predicted by our models at the 
present age ic=15 Gyr turns out to depend mainly on the 
parameters k and v of the star formation law, while being 
only slightly dependent on other model parameters. Since 
k is limited to the range 1-^2 (see § 3.2), we considered 
the extreme cases k =1, K =2 and the intermediate case 
k =1.5. To fulfill the constraint set by the current gas 
fraction, v needs to fall within a suitable range of values, 
which increase with k. For a given k, the "good" values of 
v slightly increase with the infall timescale r and with the 
mass fraction £ of short-lived stars. Here below we list our 
detailed results on the suitable ranges for v as a function of 
k, t and C- As an example, Fig. |ll] shows the evolution of 
the gas fraction in time for our Model B (k =1.5, £ =0.3, 
r =9, v =1.2, see §|9~6[). 



24 



Portinari et al. : Chemical evolution with metallicity dependent yields 



0.15 



0.1 



G(t) 



0.05 



5 10 
t(Gyr) 



15 



0.03 



a 

I 

o 
ft 



m 



0.02 



0.01 



SNII + Ib,c 



o'-i 




5 10 15 

t(Gyr) 



Fig. 11. Evolution of the normalized gas fraction G(t) for our 
model B. The vertical bar at t =15 Gyrs shows the observa- 
tional constraint on the current gas fraction 





r =5 


T =7 


r =9 


k=2.0 
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0.6-1 


0.7-1.2 


0.7-1.3 
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0.6-1 


0.7-1.2 
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9.3. The supernova rates 

The most recent empirical estimates of the Galactic SN 
rate give 14±6 (Cappellaro et al. 1997) and 21±5 (Tarn- 
man et al. 1994) SN events originating by massive stars 
- namely, Type II + Type Ib,c SNae — per 1000 yrs. 
Taking a disc radius of ^15 kpc, this roughly translates 
into (2±1)10~ 2 and (3±1)10~ 2 SNae/pc 2 /Gyr. We there- 
fore adopt (1=4)10~ 2 SNae/pc 2 /Gyr as the observational 
constraint to compare with the rate of SNae II predicted 
by our models. 

As far as SNae la are concerned, Cappel- 
laro et al. (1997) give a Galactic rate of 2±1 SN la per 
millennium, which leads to a ratio of 0.3±0.2 of Type la 
versus Type II+Ib,c SNae. Tamman et al. (1994) estimate 
that ~85% of the SNae exploding in a galaxy like our own 



Fig. 12. Evolution of the SN II+Ib,c (solid line) and SN la 
rate (dashed line) for our model B, in SNae/pc 2 /Gyr. The ver- 
tical bar at t =15 Gyr shows the observational constraints on 
the local rate of SNae II+Ib,c. The shaded area traces a factor 
0.2=0.3 of the rate of SNae II+Ib,c in time. The observational 
constraints on the relative rate of SNae la require that at the 
present age £g = 15 Gyrs the dashed line falls within the shaded 
area. 



come from massive stars, which gives a value of ~0.2 for 
the above mentioned ratio. We therefore require the num- 
ber of SNae la to be a factor 0.2=0.3 of the number of SNae 
originated by m assive stars. In our models, once Mb,i is 
fixed (see § 8.7) the efficiency with which binaries produce 
to SNae la is driven by the parameter A. For Mb,i=3 M q , 
A is constrained to be in the range 



A = 0.05 = 0.08 

independently of other model parameters. Choosing 
Mb, i=2 Mq, A would be constrained in the range 
0.03=0.05, but this alternative choice of the (Mb,i,A) has 
little effect on any practical outcome of other model pre- 
dictions. Therefore, in the following we will limit our dis- 
cussion to models with Mb : i=3 M , t4=0.05=0.08. 

As an example, Fig. |l2| shows the evolution of the lo- 
cal rate of SNae II+Ib,c and SNas la (solid and dashed 
line, respectively) as predicted by our Model B (A =0.07), 
together with the relevant observational constraints. The 
observed present relative rate of SNae la is reproduced if, 
at the present age tc=15 Gyrs, the dashed line falls within 
the shaded area, which represents a factor 0.2=0.3 of the 
SN II rate. We notice that this constraint does not sen- 
sitively depend on the assumed age of the disc, since the 
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requirement remains actually fulfilled over a long period 
(from 8 to 15 Gyrs). 

9.4- The age-metallicity relation 

The age-metallicity relation (AMR) for stars in the Solar 
Neighbourhood was first plotted by Twarog (1980). The 
data showed a large dispersion, therefore age-bins and the 
average metallicity per bin were used to trace the local 
AMR, which turned out to increase rapidly from 13 to 5 
Gyrs ago, and more slowly afterwards. Part of Twarog's 
data were re-examined by Carlberg et al. (1985), who did 
not include in their sample about fifty of the lowest metal- 
licity stars and therefore obtained a shallower slope for 
the AMR in the early phases, suggesting that the disc 
evolution started with a high initial metallicity. The bulk 
of Twarog's data was later analysed again by Meusinger 
et al. (1991), and their AMR confirmed the original result 
of a steep slope at old ages. 

A more recent study about the local AMR, based on a 
new sample of 189 nearby F and G dwarfs, was performed 
by Edvardsson et al. (1993) using high resolution spectra 
with theoretical LTE atmospheres to derive the chemi- 
cal compositions, and photometric fits with Vandenberg 
(1985) isochrones to derive ages. The resulting binned, av- 
erage AMR is in good agreement with that by Meusinger 
et al. (Fig. |p|), but the dispersion of the data about the 
average is so large that the AMR is not a tight constraint 
for chemical models. Ng & Bertelli (1997) have lately re- 
examined the Edvardsson et al. dataset, by means of mod- 
ern isochrones including new opacity tables, of new dis- 
tances from the Hipparcos dataset and giving a higher 
weight to stars with most reliable age determination. The 
resulting AMR shows a slope of ~0.07 dex/Gyr, in agree- 
ment with Edvardsson et al., but still with a large scatter. 

Additional information on the AMR comes from open 
clusters, whose age estimates are much more reliable than 
those for isolated field stars, but with the counterargument 
of the uncertainty on the correction for the radial metal- 
licity gradient and the different galactocentric distance of 
clusters. See Carraro et al. (1997) for a more comprehen- 
sive discussion of the disc AMR. 

To reproduce the large scatter in the local AMR, more 
complex models would be required, including various pos- 
sible mechanisms responsible for the scatter (orbital diffu- 
sion, non-instantaneous mixing of enriched material, self- 
propagation of star formation, local infall episodes and 
so forth). Our model applies the standard "one-zone" 
scheme, therefore it is aimed at reproducing average fea- 
tures and can only be compared to the average trend of 
the AMR, regardless of the scatter. 

The comparison with the observed binned AMR allows 
us to single out a limited range of values for the param- 
eter £. Even including the highest possible contribution 
on iron enrichment from SNae la (A— 0.08), models with 
£ <0.3 still predict too low metallicities for all the suitable 
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Fig. 13. Observational data on the local age-metallicity rela- 
tion shown together with the predictions of our Model B. 

combinations of re, v and r, and are therefore ruled out. At 
the other end, models with £ =0.5 are compatible with the 
observed AMR only if combined with the lowest allowed 
relative rate of SNae la (A=0.05); larger values for £ are 
hence ruled out, since such models would predict too high 
metallicities. Therefore, from now on we will only consider 
models with 

C = 0.3^0.5 



9. 5. The Star Formation Rate 

The present SFR in the Solar Neighbourhood is estimated 
to be ~2-10 M Q /pc 2 /Gyr (Glisten & Mezger 1982). Be- 
sides, Scalo (1986) estimates that the present SFR is 
within a factor 0.5-^1.5 of the average SFR in the past 
over the whole lifetime of the disk: 



By imposing these two observational constraints on the 
predicted SF history, we can rule out some models and 
further restrict the range of "good" combinations of the 
parameters re, v, r and £. Notably, if we assume t—5 Gyr, 
models with re=2 or 1.5 are unable to fulfill the above 
requirements on the SFR. 

Thus we further filtered the ranges of values for v listed 
here below: 
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Fig. 14. Evolution of the SFR for our model B. The dashed 
horizontal line shows the average SFR over the lifetime of 
the disc, and the vertical bar at t =15 Gyrs shows the range 
0.5-^1.5 times the average SFR. Observational constraint (see 
taxt) require that the current SFR at t — tc falls within the 
range indicated by the bar, as well as within the observational 
estimates of 2-10 M & pc" 2 /Gyr _1 
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As an example, Fig. |lj shows the evolution of the SFR 
in time for our Model B, with the related observational 
constraints. 

9.6. The G-dwarf distribution 

The G-dwarf problem, i.e. the paucity of metal-poor stars 
in the Solar Vicinity with respect to the predictions of 
the "simple closed-box model" , is usually interpreted as a 
consequence of the progressive building-up of the Galaxy 
through accretion of primordial gas, hence the failure of 
closed models in this respect (Lynden-Bell 1975; Tins- 
ley 1980). Other solutions have actually been suggested: 
prompt initial enrichment from the halo or bulge providing 
a finite initial metallicity in the ISM, higher yields due to 



Fig. 15. Observational data on the local metallicity distribu- 
tion, in relative number, of long-lived stars (dashed and dotted 
line) shown together with the predictions of our Model B (solid 
line) . 



an IMF skewed toward massive stars in the early galactic 
phases, more effective star formation in high-metallicity 
inhomogeneities of the ISM (see Pagel 1997 for a review); 
but open models with progressive inflow of primordial ma- 
terial remain the best scheme to explain the G-dwarf prob- 
lem, as well as being consistent with dynamical simula- 
tions of the formation of galactic discs and with the ob- 
servation of infall of High and Very High Velocity Clouds. 

Basing on a complete volume-limited sample of nearby 
G-dwarfs, Pagel & Patchett (1975) provided a reference 
dataset (later revised by Pagel 1989 and by Sommer- 
Larscn 1991), which was was nicely reproduced by chem- 
ical models with infall timescales of 3-4 Gyrs (e.g. Mat- 
teucci & Francois, 1989). The latest compilations of G— 
dwarf metallicities (Wyse & Gilmore 1995; Rocha-Pinto 
& Maciel 1996) show instead a narrower distribution with 
a prominent peak around [Fe/H]~-0.2, and seem there- 
fore to favour longer timescales (e.g. Chiappini et al. 1997 
and Fig. |l5|). This is also consistent with dynamical stud- 
ies, which indicate that a rather long time is needed to 
build up a galactic disc (~6 Gyrs, Burkert et al. 1992), 
especially when one takes into account that primordial or 
low-metallicity gas in the early galactic phases has a very 
low cooling efficiency (Carraro et al. 1997). 

In our models, the peak metallicity of the predicted 
distribution increases with increasing £ and A. On the 
other hand, the height of the peak mainly depends on 
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other parameters: the distribution gets narrower with in- 
creasing t, K, and v. 

All calculated models with r=5 Gyrs and t=6 Gyr 
(and with the suitable combinations of k and v as dis- 
cussed above) predict metallicity distributions of stars 
which are too broad with respect to the observed one. 
Therefore, we can rule out all these models and set the 
limit 

t > 7 

On the other hand, models with r >10 tend to predict 
too narrow distributions, therefore form now on we will 
consider only models with r =7-9 Gyrs. 

If (=0.4-0.5, whatever the iron contribution of SNae la 
(with A=0.05 to 0.08) the predicted G-dwarf distribution 
remaines peaked around [Fe/H]>0, not consistent with ob- 
servations. We rule these models out and set the limit 

C < 0.4 

In the following, we discuss models with (=0.3 as a rep- 
resentative case. Here below we list, for (=0.3, the sets of 
model parameters which are able to predict a metallicity 
distribution in agreement with the observed one, as well 
as being in agreement with all the other constraints. 
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We take as representative of the three cases k =1, 1.5 and 
2 the following models: 



model 
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0.08 
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1.5 


0.3 
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1.2 


0.07 
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We chose to use model B in Figures [El] to |l| as the rep- 
resentative model to plot versus the observational con- 
straints. In Fig. |l5|, the model well reproduces the loca- 
tion and height of the peak, while it cannot reproduce 
the higher-metallicity tail of the observed distribution. In 
fact, since our one-zone model cannot reproduce scatter 
around the average AMR, no stars are predicted to form 
with a metallicity higher than the average one reached at 
to =15 Gyrs (see also Fig |l3|) . 
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Fig. 16. Observational data on the [O/Fe] ratio vs. [Fe/H] 
as observed in the atmospheres of nearby stars. Predictions 
of model A (dotted line), B (solid line) and C (short-dashed 
line) are shown for comparison. ^'The data from Tomkin et 
al. (1992) are taken as revised by Carretta (1995) 



9.7. Abundance ratios 

Figures |l6| to ^ display observational data on the abun- 
dance ratios of various elements with respect to iron, to- 
gether with the predictions of our models A, B and C. The 
data refer to spectroscopic abundances in the atmospheres 
of nearby stars; detailed data sources are indicated in the 
legends of the figures. 

Predictions for the a-elements ( 16 0, 24 Mg, 28 Si, 32 S, 
40 Ca) seem in most cases to well reproduce the trend of 
the data. A remarkable exception is the [Mg/Fe] ratio. 
Our yields for the SN explosion are derived from those of 
WW95, and therefore display the same problem of 24 Mg 
underproduction (Timmes et al. 1995). As discussed by 
Thomas et al. (1997), uncertainties in SN calculations still 
allow for rather different yields among different authors, 
especially in the case of 24 Mg production. For instance, 
24 Mg yields by Thielemann et al. (1996) are sensitively 
higher and better suited at explaining the [Mg/Fe] ratio 
observed in low-metallicity nearby stars. See, in this re- 
spect, the lower panel of Fig. |23|, where we compare the 
24 Mg ejecta from the CO-cores by WW95 with those 
by Thielemann et al. (1996) (these latter obtained in a 
similar fashion as for WW95, see App. A. As in Timmes 
et al. (1995), there is a 32 S underproduction which could 
not be solved by adopting Thielemann et al. 's yields, 
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Fig. 17. Observational data on the [Mg/Fe] ratio vs. [Fe/H]. 
Model predictions are plotted as in Fig. 16. (1) The data from 
Zhao & Magain (1990) are taken as revised by Carretta (1995) 
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Fig. 18. Observational data on the [Si/Fe] ratio vs. [Fe/H]. 
Model predictions are plotted as in Fig. 16 



which in the case of sulfur are even lower than those by 
WW95. 

Our model predicts an overproduction of carbon with re- 
spect to the observational data. The mismatch with the 
data starts at low metallicities, indicating that the con- 
tribution of short-lived, massive stars to C-production is 
overestimated. Timmes et al. (1995) do not display a sim- 
ilar problem when using directly WW95 SN ejecta. At 
least two reasons might be tracked for this. (1) Stellar 
models with overshooting give a lower limiting mass for 
SN explosion (M up ~6 Mq rather than 8 M©) and there- 
fore C-producing massive stars have a larger global weight 
in the IMF is taken into account. (2) Our link between 
mass-losing pre-SN structures and SN models requires 
some assumptions, especially relevant to the edge of the 
He-burning region, where carbon is produced; this might 
make 12 C yields particularly uncertain in our grid (see 
App. B). With respect to the data, our carbon overpro- 
duction in massive stars amounts to a factor of ~2, which 
is still acceptable when taking into account the uncertain- 
ties in SNae II yields. Again, if for instance we compare 
the results by WW95 and by Thielemann et al. (1996) for 
a given Mco, the latter give 12 C yields up to a factor of 
2 lower (see Fig. 23, upper panel). 

Unlike for the a elements, the predicted [N/Fe] ra- 
tio at low metallicities depends sensitively on the past 
SF history. This effect is mainly due to the contribution 
of primary 14 N from intermediate mass stars, which in 
the range 4-5 M© undergo envelope CNO-burning during 
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Fig. 19. Observational data on the [S/Fe] ratio vs. 
Model predictions are plotted as in Fig. 16 
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Fig. 20. Observational data on the [Ca/Fe] ratio vs. [Fe/H] 
Model predictions are plotted as in Fig. 16 
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Fig. 22. Observational data on the [N/Fe] ratio vs. [Fe/H]. 
Model predictions are plotted as in Fig. 16. The long-dashed 
line shows separately the contribution of the sole secondary 
component for Model B (see text) 
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Fig. 21. Observational data on the [C/Fe] ratio vs. [Fe/H]. 
Model predictions are plotted as in Fig. 16. (1) The data from 
Tomkin et al. (1992) are taken as revised by Carretta (1995) 



the TP-AGB phase (Marigo et al. 1997). These stars re- 
lease their contribution of primary 1 N on timescales of 
^1 — 2 10 8 yrs, independently of the initial metallicity. In 
Fig. ^2] we display the relevant data versus the the pre- 
dictions of our models A, B, C. As k decreases (from C 
to A), the SF gets less efficient in the early phases, and 
the metallicity [Fe/H] evolves more slowly. Therefore, the 
[N/Fe] ratio increases more rapidly in the early phases if 
k is lower. 

Abundance ratios of a elements do not show a similar 
dependence on the SF efficiency in the early phases, since 
in the first Gyr or so iron and a elements are produced by 
the same source (massive stars), and their relative abun- 
dance is determined by the adopted yields and IMF for 
massive stars, irrespectively of the SF history. On the con- 
trary, in the case of 14 N the effect is related to the compe- 
tition between different sources for the two elements: if k is 
higher, more gas is recycled through successive generations 
of massive stars in the first ~ 10 8 yrs, and intermediate 
mass stars contribute to the rise of the [N/Fe] ratio when 
the metallicity [Fe/H] is already relatively high. In Fig. 
we also display the contribution of secondary 14 N alone 
(long dashed line); unlikely the primary component, the 
[N/Fe] ratio due to the bare secondary component is basi- 
cally independent of the SF history, therefore we just show 
Model B for the secondary component. Evidently, primary 
14 N from intermediate mass stars is a major component 
in the overall 14 N production, and must be properly taken 
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Fig. 23. Ejecta of 12 C (upper panel) and 24 Mg (lower panel) of 
the CO-core vs. Mco as deduced from Thielemann et al. 1996 
(open circles) compared with those from WW95 (case Z =Zq, 
filled circles) 

into account. Still, even including this component, models 
seem to underproduce 14 N with respect to 56 Fe, especially 
in the low-metallicity range. A classical suggestion to over- 
come this mismatch is to assume some primary 14 N contri- 
bution by massive stars (Matteucci 1986). This might be 
possible in very low-metallicity massive stars if some pecu- 
liar mixing event is allowed to occur (Timmes et al. 1995). 
Fig. ^ might also suggest, as an alternative, that a slow 
SF history (n —\) combined with a more efficient produc- 
tion of primary 14 N in intermediate mass stars might solve 
the problem. A more extreme alternative might be an ini- 
tial IMF skewed to massive stars, which quickly provide 
the ISM with a minimum amount of CNO isotopes so that 
secondary production of 14 N becomes effective already at 
very low metallicities. But this alternative should need to 
be tested also with respect to the other abundance ratios 
and the other model constraints. 

10. Summary and conclusions 

We have presented a detailed study of the ejecta and chem- 
ical yields from stars of different mass and initial chemical 
composition, taking into account modern results of hydro- 
static and explosive nucleosynthesis. The analysis is made 
for stars that avoid the AGB phase, i.e. from 6 M up 
1000 M Q . 

In this range of mass wc distinguish three main groups: 
(i) the 6 to 9 M stars, which are not significantly affected 



by mass loss by stellar wind and terminate their evolution 
by the electron capture instability and consequent super- 
nova explosion; (ii) the 9 to 120 M stars, which suffer 
from mass loss by stellar wind starting from the zero age 
main sequence and terminate as iron-core collapse or PC 
supernovae; and finally (iii) the rather hypothetical range 
of very massive objects (120 to 1000 M ), which are char- 
acterized by pulsational instability and violent mass loss 
during the core H-burning phase and may terminate the 
evolution either as iron-core collapse or "PC" supernovaeor 
black holes, depending on the size of their CO-core 

We start from the sets of hydrostatic stellar models 
by the Padua group that are calculated in presence of 
mild convective overshoot. This explains why the upper 
mass limit for the occurrence of the AGB phase is set at 
about 5 Mq depending on the initial chemical composi- 
tion. All these evolutionary sequences are carried out from 
the ignition of the core H-burning till the start of central 
C-burning, taking into account mass loss by stellar wind 
wherever appropriate. These sets of models span a wide 
range of initial masses and chemical composition and are 
fully homogeneous as far as the input physics, accuracy 
and numerical methodology are concerned. So they are 
the ideal tool to construct equally homogeneous sets of 
ejecta and stellar yields. 

From the models we derive the basic relationships be- 
tween the initial and final (namely, at the onset of C- 
burning in the core) mass of the stars, and between the 
initial mass and the final mass of the helium and carbon- 
oxygen core. The subsequent evolutionary stages are so 
short-lived that the above relationships will not change at 
all. If mass loss occurs by stellar wind during the core H- 
and He-burning stages, the amounts of mass in the various 
species lost in the wind are calculated. 

Adopting the CO core mass as the driving parameter 
we link the hydrostatic structure at early core C-burning 
to the final explosive stage. Using literature data on su- 
pernova models, we calculate the amount of mass ejected 
by each CO core at the time of explosion for each species 
under consideration, to which the contribution from the 
overlying layers still existing in the star at this stage, and 
the contribution from stellar wind in the earlier stages are 
added. 

The results are summarized in a series of tables, and 
all details relative to construction of the ejecta and yields 
for massive and very massive stars are given in the main 
text and Appendices A and B. 

The corresponding ejecta and yields for low and inter- 
mediate mass stars (those passing though the AGB phase) 
can be found in Marigo et al. (1996, 1997). They are based 
on the same stellar models in usage here up to the end 
of the E-AGB, and are derived from an original semi- 
analytical method suitably tailored to follow the TP- AGB 
phase in presence of envelope burning in the most massive 
stars of this group. So the ejecta and yields for massive and 
very massive objects presented here together with those 
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by Marigo et al. (1996, 1997) constitute a dataset of un- 
precedented homogeneity. 

Our results about stellar nucleosynthesis and ejecta 
have been applied to study the chemical enrichment of 
the Solar Neighbourhood by means of the infall model of 
Chiosi (1980), in which the matrix formalism developed 
by Talbot & Arnett (1973) to calculate the fraction of a 
star of initial mass M that is ejected in form of chemi- 
cal species i has been incorporated. This matrix formal- 
ism is particularly useful to understand what fraction of a 
star originally in form of species j is eventually ejected as 
species i. The numerical method to solve the set of model 
equations governing the chemical evolution of the ISM is 
the same as in Talbot & Arnett (1971), however adapted 
to the infall scheme. 

The chemical model contains a number of parameters, 
some of which are directly fixed by the observations or 
estimated independently (e.g. the distance of the Sun from 
the Galactic center, the local surface mass density and the 
age of the galactic disk) . Other parameters are let vary so 
that some basic observational constraints are met. All our 
results refer to the Scalo IMF with fixed slope and upper 
mass limit M„=100 M Q . 

Constraints of the models are (i) the current gas frac- 
tion; (ii) the rate of type I and type II SNae; (iii) the age 
metallicity relation; (iv) the past and present estimated 
star formation rate; (v) the distribution of long-lived stars 
in metallicity (G-dwarf problem) . 

Free parameters of the models are: (i) the efficiency 
and exponent of the SFR {v and k, respectively); (ii) the 
infall time-scale r; (iii) the mass fraction ( of the IMF 
in stars with M >1 M Q ; (iv) the amplitude factor A and 
minimum mass of the binary Mb,i for Type la SNae. 

A careful comparison of model results with the obser- 
vational constraints allows us to disentangle the effects 
of different parameters on the model predictions and put 
successive limits on them. We are finally left with the fol- 
lowing "good ranges" for the various parameters: 1< k <2, 
0.3< v <3 (in suitable combination with k), 0.3< £ <0.4, 
7< r <9 Gyrs, 0.05< A <0.08 (for M B ,i='i M ). 

Predictions for most elemental ratios as a function of 
[Fe/H] seem to reproduce the trend of the data, but for 
[Mg/Fe] and to a less extent for [C/Fc] and [N/Fc]. Pos- 
sible reasons for this disagreement are briefly discussed. 
In the case of 24 Mg and 12 C, they might be reconduced 
to uncertainties in the ejecta of SNae. In the case of 14 N, 
the predicted [N/Fe] is crucially dependent on the past 
SF history due to the contribution of primary 14 N from 
intermediate mass stars. Possible ways out to the low pre- 
dicted 14 N abundance are briefly recalled: either primary 
14 N from massive stars of low metallicity (Matteucci 1986) 
or an IMF more skewed toward massive stars than pre- 
dicted by a Salpeter-like law during the very early stages 
of galaxy formation and evolution (cf. Chiosi et al. 1997). 
Both hypotheses need to be thoroughly investigated. The 
problem is open. 



Finally, we like to remark our effort to derive metallic- 
ity and hence time dependent stellar yields, since limiting 
to consider constant yields is not compatible with our cur- 
rent understanding of stellar evolution. But we also need 
to underline that for some elements many uncertainties 
still exist on the nucleosynthetic yields of stars in the dif- 
ferent mass ranges. These uncertainties display already 
when considering model predictions for the Solar Neigh- 
bourhood, and should be borne in mind when applying 
chemical models to external galaxies, where detailed tests 
of the predictions are more difficult. 
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Appendix A: the link with SN models 

Here we describe in detail the method we adopt to establish the 
relations between Mco and the amount of different elements 
expelled by the core, so that we can "complete" the total ejecta 
of our models with the outcome of the final iron-core collapse 
SN. To this aim, we need to deduce the contribution of the sole 
CO-core Ef° to the total ejecta of the SN models by WW95. 
In the following, we adopt this notation: 

M total mass of the model by WW95 

Mco mass of the CO-core of the model 

Z metallicity of the model 

Ei total ejected amount of species i, taken from 
WW95 tables 

E" ew newly synthesized and ejected amount of species 
i 

Xf initial abundance of species i in the model 
Ef xt amount of species i contained in the layers exter- 
nal to Mco 

Ef° amount of species i ejected by the CO-core, i.e. 

contribution of the CO-core to the global ejecta 
of i 

M env mass of the region unaffected by CNO burning 
(roughly corresponding to the mass of the enve- 
lope, apart from dredge-up episodes) 

First of all, for each tabulated total mass we need to de- 
duce the relevant value of Mco- We assume that the expelled 
amounts of H and He originate solely in the layers overlying 
Mco, that is to say (1) in the H-rich envelope, possibly en- 
riched in 4 He and in 14 N by convective dredge-up episodes, and 
(2) in the fraction of the He-core that was not involved in the 
He — > CO burning, Mue- Mco- (With these assumptions we 
neglect a possible contribution in He from an a rich freeze- 
out, which anyway should contribute only a small fraction of 
the total He ejecta.) The global metallicity in the layers exter- 
nal to Mco is assumed to correspond to the initial one, since 
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in these layers no nuclear burning takes place but H-burning 
and the CNO cycle, though altering the relative abundances of 
CNO isotopes, leaves the overall metallicity substantially un- 
changed. Therefore, on the base of the expelled amount of H 
and He, and correcting for the initial fraction of metals, we can 
derive the mass of the layers over the CO-core; then we derive 
Mco as: 



M CO = M 



E H + E H e 
1 - Z 



(Al) 



In the 15 and 25 M Q , Z=Zq cases (Tab. S in WW95), we get 
Mco = 2.53 and 6.59 Mq respectively; these values look rea- 
sonable when compared to Figs. 10a, b of WW95, that illustrate 
the chemical profiles of the corresponding pre-SN structures. 
This gives us confidence in the method we adopt (but see § [l(]) . 
Once we get Mco we consider the production sites of the vari- 
ous elements, so that we can distinguish the ejecta originating 
in the CO-core from the contribution of the overlying layers. 

14 N — 13 C are produced in the CNO cycle and destroyed by 
a capture during He-burning. No 14 N and 13 C are left in- 
side the CO-core, and they are not produced by explosive 
nucleosynthesis (see also Tab. 7 in WW95, comparing the 
chemical composition of a 25 Mq model before and after 
explosion). These species can come only from layers out- 
side Mco , which doesn't contribute to their ejecta; see also 
Fig. 13a,b of WW95. 

Eaix = 



J C13 
pCO 







E N 

12 C — 16 O are synthesized during He-burning inside Mco - In 
order to remove the contribution of the outer layers, we 
need to estimate carbon and oxygen abundances outside 
Mco- In layers experiencing the CNO cycle, their original 
abundances are altered in favour of 14 N and 13 C. Conse- 
quently, we estimate first the amount of newly synthesized 
14 N and 13 C. Since these are located only outside Mco, 
by subtracting the initial amounts within M-Mco to their 
total ejecta we directly get the newly synthesized amounts: 

E c Tz = Ecu - X C13 (M - Mco) 
Eff w — Em — X%(M — Mco) 

13 C is synthesized in the CNO cycle starting from 12 C, 
while 14 N is synthesized starting from 12 C and 16 0. As- 



suming that the new 14 N originates from 12 C and 10 O pro- 



le, 



portionally to their initial abundances, the amount of 12 C 
and 16 O lying outside Mco can be expressed as: 



E e c xt = X C (M - Mco) 
E Q xt = X° (M - Mco) ~ E 



j->new 



ti N 



Xr 



X°c + X" 



X° c + X 

Here, the first term is the amount of 12 C and 16 O present 
outside Mco since the beginning, from which we remove 
the amount that has been converted into 14 N and 13 C. The 
contribution of the CO-core to carbon and oxygen ejecta 
is: 

Ec° = Ec 



Eo° = Eo 



FjO 



15 N is quickly destroyed by the CNO cycle, which reduces 
its abundance by an order of magnitude. The SN explo- 
sion can produce 15 N by neutrino nucleosynthesis on 16 O 
( 16 0(v x ,v x p) 15 ~N). We assume that the abundance of 15 N 
is negligible where the CNO cycle takes place, unaltered 
with respect to the initial one elsewhere. We estimate the 
mass not involved in H-burning as: 

M - — 

Starting with the total ejecta of 15 N, we subtract the con- 
tribution of the "envelope" (mass unaffected by the CNO 
cycle) with the initial abundance of 15 N. Whatever is left 
comes from neutrino nucleosynthesis on 16 O, both on the 
16 O in the CO-core and on the 16 O in the outer layers; we 
further scale the resulting amount of 5 N with respect to 
the fraction of 16 O within Mco - The resulting contribution 
of Mco to the ejecta of 15 N is expressed as: 



-C-jvis 



(En15 — Menv X 



e; 



N15) 



O 



Eo 



17 abundance within Mco is negligible because 17 is de- 
stroyed by a capture during He-burning — see also 
Figs. 13a,b of WW95. 

Eo?7 = 

18 O is destroyed in the CNO cycle and later produced, but 
also rapidly destroyed, during He-burning; no ls O is left 
inside the CO-core. 



20 Ne— 24 Mg— 28 Si— 32 S 



P co 



40 Ca are produced during He- 
and C-burning; we assume that out of the CO-core their 
abundances are unaltered with respect to the initial ones. 
To get the contribution of the CO-core, we simply remove 
from the total ejecta the contribution of the overlying lay- 
ers, where the initial abundances hold: 

(M-M CO )X° Ne 

tV : = E Mg - (M - Mco) X° Mg 

Esi 1 = E,si — (M — Mco) X ^ 

E c s ° = E S - (M - Mco) X° s 

E§° =E G a-(M-Mco)X Ca 

50 Fe is produced in the very last hydrostatic burning stage (Si- 
burning), but most of the ejected amount of iron originates 
in the radioactive decay of 56 Ni. The tables from WW95 
list the ejecta at a time immediately after explosion, when 
the decay hasn't taken place yet; since all the released 56 Ni 
should later decay in 56 Fe, we add the ejecta of 56 Ni to 
those of 56 Fe: 

-(M-M co )X° Fe ] + E m 



? co 

^CO 



pCO 



[E f 



In the case of 30, 35 and 40 Mq WW95 calculate several mod- 
els (A, B, C) differing in the energy of the ejected material at 
infinity. This results mainly in different ejecta of 56 Fe, since 
this element is sensitive to the location of the mass cut, which 
changes with the explosion energy. In our link we referred to 
the A models of WW95; the effect of different assumptions 
about the explosion energy can be included in the uncertainty 
of a factor of two upon the amount of ejected iron, mentioned 
by WW95. 
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Appendix B: some warnings about the link 

Here we discuss possible drawbacks or incoherences of the 
method we adopted to link our stellar models with WW95 
SN models. 

(1) We have implicitly assumed that explosive nucleosyn- 
thesis only involves the CO-core leaving the overlying layers 
unaffected, as they are added from our pre-SN models. This 
sounds a reasonable approximation, but we neglect a possi- 
ble production of 15 N from neutrino nucleosynthesis over the 
16 O lying out of Mco and a possible contribution of an a rich 
freeze-out in the ejecta of 4 He. 

(2) One of the basic parameters in the adopted SN mod- 
els is the kinetic energy of the ejecta at infinity (WW95): the 
Shockwave generating the explosion is simulated by means of 
a piston, regulated so that the final energy assumes a typical 
value, ~ 1.2 x 10 51 erg. Outside Mco, the layers in our pre- 
SN models have a different mass and structure from those of 
WW95, both because we include mass loss and because of the 
different physical treatment of details of stellar evolution (es- 
pecially the inclusion of convective overshooting). Therefore, 
for a given Mco, the kinetic energy of our ejecta at infinity 
can be different from that imposed in WW95 models. This 
affects mainly the remnant mass, which depends on the explo- 
sion energy, on the pre-SN structure, on the stellar mass and 
on metallicity (WW95). Most of all, the extent of a possible 
reverse shock as found in WW95 models may not be directly 
transferred to our models, because it depends on the density 
structure of the layers over the core. A reverse shock can in- 
duce a fall back of material toward the collapsed core, which 
increases the mass of the final remnant. Since we do not know 
how the remnant mass in WW95 models depends on this physi- 
cal effect in detail, we can't correct for it. The choice of adopt- 
ing the same remnant mass for the same Mco as in WW95 
models can influence our ejecta of 28 Si, 32 S and 56 Fe, that arc 
the elements produced next to the mass cut of the SN. 

(3) An inconsistency between the two sets of models is the 
different cross-section adopted for the fundamental reaction 
12 C(a, 7 ) 16 0. In our tracks the adopted value is from Caugh- 
lan & Fowler (1988), while WW95 adopt 1.7 times that value. 
Both values are within current uncertainties, but they result in 
a sensitively different composition of the CO-core: in our mod- 
els the resulting core consists of a 30% of 12 C and a 70% of 
16 0, while in WW95 models the carbon abundance in the core 
is as low as ~10% (see Figs. 5ab, lOab, 12ab, 14ab and 16ab of 
WW95). In a way, by "replacing" our CO-cores with the cores 
of WW95, we adopt their value of the cross-section, since that 
reaction is important only from He-burning on, i.e. inside the 
CO-core. However, in cases where mass loss is so efficient as 
to reveal He-burning processed material on the surface and 
disperse it in the wind, the carbon/oxygen abundance ratio in 
the expelled material is due to a cross-section value which is 
different from that used inside the core; this gives some incon- 
sistency in case of extreme mass- loss. 

(4) Most of all, our method neglects the finite extension of 
the He-burning shell: when we separate the CO-core contribu- 
tion from the contribution of the overlying layers, we somehow 
assume a sudden jump in the chemical composition between 
the CO-core and the outer layers, as if He-burning defined 
Mco as a sharp edge. Actually, there is a gradient of chemical 
profile between the regions where helium is still abundant and 



those where it has been completely processed. In cases where 
the He-burning shell is not thin and/or a shallow gradient 
of chemical composition has been established, our assumption 
about a clean separation between the CO-core and the rest 
of the star may lead to inconsistency. This would mainly af- 
fect the early products of He-burning, such as carbon, whose 
abundance can change sensitively even in the outer regions of 
He-burning, where the process is only partial. This brings some 
scatter in the C-production vs. Mco relation, since at the time 
of explosion most of the carbon is located at the edge of the 
CO-core or even outside it, and its distribution is sensitive 
to past history especially in stars which suffered from efficient 
mass loss (M92). Therefore, the linking method can be less reli- 
able in the case of carbon, and this might explain our problems 
in reproducing the observed [C/Fe] ratio (see § 9.7). To check 
this effect, we apply the method outlined in App. A to derive 
the ejecta of the CO-cores of the mass-losing helium stars cal- 
culated by Woosley et al. (1995). These stars build up a shal- 
low gradient of carbon and helium abundance at the outer edge 
of the CO-core, and for these stars the assumption of a sharp- 
edged core is not properly correct. By applying Eq. (Al), we 
tend to overestimate the size of the CO-core for these stars, 
as shown in Tab. |l5| where the estimated core masses are com- 
pared with the real ones as tabulated in Woosley et al. (1995, 
Tab. 6). In Fig. 24 we plot the ejecta versus Mco for these he- 
lium stars and we compare these relations to those obtained for 
WW95 models of Z=Z Q (App. A and Fig. §). While carbon 
evidently behaves in quite a different way in the two cases, the 
ejecta of other elements and the remnant mass seem to follow 
a better defined relation with Mco- 



Table 15. I s * column: model masses of mass-losing he- 
lium stars of Woosley et al. 1995; 2 nd column: corresponding 
CO-core masses as listed in Tab. 6 of Woosley et al. 1995; 3 rd 
column: CO-core masses as deduced according to the outlined 
method (see text) 



M 


Mco(l) 


M CO (2) 


4 


1.53 


1.60 


5 


1.87 


2.05 


7A 


2.30 


2.73 


7B 


2.30 


2.80 


7C 


2.30 


2.80 


10A 


2.50 


3.26 


10S 


2.50 


3.26 


20 


2.53 


3.31 



Appendix C: calculating the Qij matrix 

Here we explain the definition and the meaning of the non- 
zero elements of the Qij matrix (Tab. p^[ ) and show how the 
fundamental quantities entering the matrix are calculated. 

1 H For any stellar mass, all the hydrogen which remains un- 
burnt in the interiors is eventually ejected. The only non- 
zero matrix element for *H is Qi,i, corresponding to the 




Fig. 24. Remnant mass and ejecta of the CO-core vs. Mcoas deduced for mass-losing helium stars of Woosley et al. 1995 
(open circles) compared with the analogous relations deduced for the case Z=Z & of WW95 (filled circles). 
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initial 1 H which hasn't been processed. As §4 is the mass 
fraction where 1 H is turned to 4 He, 

Qi.i = 1 — <74 

2 He Helium is synthesized on 1 H within q± and later turned to 
12 C, 16 and heavier species within qc- All the surviving 
4 He (at least for stars of some interest, M >0.5 Mq) is 
eventually ejected. The term corresponding to the newly 
synthesised and ejected 4 He is 

Q2,i = (74 - qc 

while the term corresponding to the original unprocessed 
and ejected 4 He is 

Q2.2 = 1 - qc 

12 C Carbon is turned to 13 C and later to 14 N within qci3s, by 
the CNO cycle; later, it is synthesised on 4 He within qc- 
The original 12 C which is ejected is represented by 

Q3,3 = 1 — qci3s 

Also the newly synthesized 12 C within wc = qc — d is not 
locked in the remnant and can be ejected: 

Qi,i = <?3,2 = Xcwc 

The matrix element Qs,i is non-zero because the 3a chain 
turns to 12 C both the original 4 He and the 4 He which had 
been previously synthesized on the original 1 H. 
13 C The ejected original 13 C is the one outside qNs, where it 
is turned to 14 N by the CNO cycle: 

<?4,4 = 1 — qNs 

Some new 13 C can also be synthesized on 12 C by the CN 
cycle within qci3s, and the corresponding ejected fraction 
is represented by: 

Qi,3 = qcvia — qNs 

In intermediate mass stars some primary 13 C can also be 
produced, due to the combined effect of III dredge-up and 
CNO envelope burning, so for these stars further non-zero 
matrix elements can be: 

Q4,l = £?4,2 = XC13WC 

14 N Only the original 14 N outside the He-burning region can 
be eventually expelled, otherwise it is turned to neutron- 
rich isotopes: 

Qs.s = 1 - qc 

Secondary 14 N is synthesized on 12 C, ie O and 13 C during 
the CNO cycle; the fraction which remains outside qc is 
eventually ejected for any stellar mass: 

Qb,3 = Q5,4 = Q5,6 = QNs — qc 

We point out that, by taking Qs,3 = Q5.4, = Qs,6, we im- 
plicitly assume that in the CNO cycle secondary 14 N comes 
from original 12 C, le O and 13 C in the same proportions as 
their initial abundances. 

In intermediate mass stars some primary nitrogen is pro- 
duced by CNO envelope burning over dredged-up 12 C, so 
for these stars we also have: 

Qs.i = Qs,2 = Xnwc 
16 The ejected original oxygen is that outside <jjv s , where it 
is turned to 14 N by the ON cycle: 

Qe,6 = 1 — qNs 

New oxygen is produced within qc; the fraction which isn't 
locked in the remnant d is expelled: 

Qe,i = Qe,2 = Xowc 



nr We indicate with nr the body of neutron rich isotopes 
synthesized on 14 N during He-burning ( 18 , 22 Ne, 25 Mg) 
within qc- All the original and newly synthesized nr iso- 
topes are eventually ejected; the corresponding matrix ele- 
ments are: 

<?7,7 = l-d 

Q"7,3 = QlA = <?7,5 = Ql,6 = WC 

20 Ne — 24 Mg— 28 Si— 32 S— 40 Ca — 56 Fe All these species are 
produced in the He-burning stage or in later stages, within 
qc- The fraction of the original amount of each of these 
species which is eventually ejected is represented by: 

Qi,i = 1 - d i = 8 -M3 

while the ejected fraction of the newly synthesized amount 



Qi,i = Q 



i,2 = Xi W C 



8^13 



The relevant quantities entering the Qij matrix (§ 3.4) are 



calculated on the base of the stellar ejecta derived from the 
tracks of the Padua library and discussed in this paper. We 



impose, species after species, Eq. 
the proper abundance set {Xj , j - 
of stellar ejecta are preserved. 



so that by re-inserting 
..13} the original values 



d 



d = 



(CI) 
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The remnant mass fraction d is, by definition: 

M r 

~M 

By definition, 54 is the mass fraction involved in the X H^ 
4 He burning. If we simply considered the mass fraction of 
the He-core, 54 = Mhs/M, we would neglect that the core 
may have receded or varied in size in the course of evolu- 
tion and that some material has been dredged-up in the 
envelope and/or lost in stellar winds. On the contrary, by 
applying Eq. ( |^ ) to the case of hydrogen (i=l): 

Eh 
~W 

we get an "effective" 54 which automatically takes into ac- 
count the whole mass where hydrogen has been burnt, both 
the mass lying in the core and the mass dredged-up to the 
surface. So we assume: 
E H 



(l-g 4 )X 



?4 



XM 



(C2) 



qc By definition, qc is the mass fraction where He is turned 
to 12 C, 16 and heavier species. Again, if we simply took 
the mass fraction of the CO-core, qc = Mco /M, we would 
neglect the material dredge-up in the envelope and/or 
ejected with mass loss. In the case of low and intermediate 
mass stars, we would lose completely the contribution of 
the III dredge-up, since we'd find qc = d and wc = 0. On 
the contrary, by applying Eq. (^) to the case of helium 
(i=2), we get an "effective" qc which automatically takes 
into account the role of dredge-up episodes: 

Ehb 



( qi -qc)X + (l-qc)Y 



qc 



q 4 X + Y 



Eh. 
M 



X + Y 



(C3) 



where 54 is known from Eq. (C2). Now we can also deter- 
mine: 

w c = qc ~ d (C4) 
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gjvs We obtain the mass fraction qNs where secondary nitro- 
gen is produced by imposing Eq. (^3|) to the case of sec- 
ondary 14 N. This corresponds to all the ejected 14 N in 
the case of massive stars, while in low and intermediate 
mass stars we keep the ejecta of primary and secondary 
14 N apart, since they involve distinct matrix elements. 

{qNs ~ qc)(X c + X C i3 + Xo) + (1 - qc)X N = ^ 



Ens 



(Xc + Xcvi + Xo)M 
Xn 



(1-flo) 



Xo + Xoia + Xo 



+ qc 



(C5) 



where qc is known from Eq. (S3). 
qci3s We obtain the mass fraction qciss where the CN cycle 
turns 12 C to 13 C in the same way as qNs, i.e. by imposing 



Eq. ( |23j ) to the case of secondary C: 
(qci3s — qN 3 ) Xc + (1 — qNs) Xci3 



Eci 



M 



qci3s — 



Eci 



It -,Xci3 . 

- (1 - qNs)— — h qNs 



(C6) 



X C M v X c 
where qNs is known from Eq. (C5). 
Xjv The abundance \n of the 14 N synthesized within wc by 
envelope burning in intermediate mass stars is defined by 
Eq. (E3), applied to primary 14 N: 
En p 



(Xnw c )(X + Y) = 
En p 



M 



XN 



(X + Y)Mw c 
where wc is known from Eq. (\ 
Xci3 Similarly, by imposing Eq. (E 

Ec\3 V 

XC13 - 



(C7) 

53). 

31) to primary 13 C we get: 
(C8) 



(X + Y)Mw C 

Xc The abundance xc of new 12 C synthesized within wc is 



obtained by applying Eq. ( 

(Xcw c )(X + Y) + (1 - q ci3 s) X c 
Ec 



Xc = 



1 - qci3 



to carbon («=3) 

= E° 

M 
X c 



(X + Y)Mwc wc 
where qci3s is known from Eq. 
Xo Similarly, by applying Eq. @ 
Eo 1 - qNs 



Xo 



X + Y 

to 16 O (i=6) we get 
Xo 



(C9) 



(CIO) 



( Xl wc)(X + Y) + {l-d)X 1 = 



(X + Y)Mw c w c X + Y 
where qjvs is known from Eq. (C5). 
XAfe — XMg — XSi — Xs — XFe The abundances of the newly 
synthesized component of these elements within wc is ob- 
tained by imposing Eq. J23I ) to the cases i=8-^13. 

Ei 
M 

Ek 1-d Xj 

Xl (X + Y)Mw c wc X + Y 1 ' 

We can now derive a Qij matrix for each star whose ejecta 
have been determined; then, using Eq. dill) , we get the "resti- 
tution fractions" Rah entering Eq. (7). We have derived stel- 
lar ejecta and Qij matrices for 5 different sets of metallicities; 
therefore in Eq. (^l|) we have an explicit dependence on the 
initial composition of the star through the Xj's, and also an 
indirect dependence through Qij = Qij(M, Z). That's why the 



restitution fractions Rkh's in Eq. (7) are to be evaluated as a 
function of the birth-time of the star, to take the effect of initial 
composition into account: 



i(t - t m ) = Qii(M, Z(t - r M ))Xj(t - t m ) 



It is worth here underlining the advantages of this u Qij matrix" 
formalism. It was originally introduced in order to compensate 
for the lack of stellar model of different metallicities, by as- 
suming that the influence of the initial chemical composition 
over the -Rmi's was included in the linear dependence on the 
Xj's, while the dependence of the Qij's on Z was negligible 
(Talbot & Arnett 1973). Now that stellar models and corre- 
sponding ejecta for different metallicities are available, using 
the Qij matrix allows to take into account also possible differ- 
ences of chemical composition within a given Z. Indeed, stellar 
models with different metallicities generally assume solar rel- 
ative abundances of the various species within a given Z; but 
abundance ratios are not constant in the course of galactic evo- 
lution, nor they are in the evolution of a chemical model. The 
Qij matrix links any ejected species to all its different nucle- 
osynthetic sources, allowing the model to scale the ejecta with 
respect to the detailed initial composition of the star through 
the X/s. 
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